International Journal of Modern Physics D 
© World Scientific Publishing Company 



m 

o . 

, On the structure of the burst and afterglow of Gamma-Ray Bursts I: 

' the radial approximation 

D 

^ ■ REMO RUFFINI, CARLO LUCIANO BIANCO and SHE-SHENG XUE 

ICRA — International Center for Relativistic Astrophysics and Dipartimento di Fisica, 
Universita di Roma "La Sapienza", Piazzale Aldo Mora 5, 1-00185 Roma, Raly. 

PASCAL CHARDONNET 

ICRA — International Center for Relativistic Astrophysics and Universite de Savoie, 
LAPTH - LAPP, BP no, F74941 Annecy-le-Vieux Cedex, France. 

O ! FEDERICO FRASCHETTI 

ICRA — International Center for Relativistic Astrophysics and Universita di Trento, 
Via Sommarive 14, 1-38050 Povo (Trento), Italy. 



o 



Oh' 



o 



We have recently proposed three paradigms for the theoretical interpretation of gamma- 
ray bursts (GRBs). (1) The relative space-time transformation (RSTT) paradigm em- 
phasizes how the knowledge of the entire world-line of the source from the moment 
^ ' of gravitational collapse is a necessary condition in order to interpret GRB datalD (2) 

The interpretation of the burst structure (IBS) paradigm differentiates in all GRBs 
between an injector phase and a beam-target phase. ^ (3) The GRB-supernova time 
. , - sequence (GSTS) paradigm introduces the concept of induced supernova explosion in 

' the supernovae-GRB association0 In the introduction the RSTT and IBS paradigms 

?H are enunciated and illustrated using our theory based on the vacuum polarization pro- 

cess occurring around an electromagnetic black hole (EMBH theory). The results are 
summarized using figures, diagrams and a complete table with the space-time grid, the 
fundamental parameters and the corresponding values of the Lorentz gamma factor for 
GRB 991216 used as a prototype. In the following sections the detailed treatment of the 
EMBH theory needed to understand the results of the three above letters is presented. 
We start from the considerations on the dyadosphere formation. We then review the 
basic hydrodynamic and rate equations, the equations leading to the relative space- time 
transformations as well as the adopted numerical integration techniques. We then illus- 
trate the five fundamental eras of the EMBH theory; the self acceleration of the e+e~ 
pair-electromagnetic plasma (PEM pulse), its interaction with the baryonic remnant of 
the progenitor star, the further self acceleration of the e+e~ pair-electromagnetic ra- 
diation and baryon plasma (PEMB pulse). We then study the approach of the PEMB 
pulse to transparency, the emission of the proper GRB (P-GRB) and its relation to the 
"short GRBs" . Particular attention is given to the free parameters of the theory and to 
the values of the thermodynamical quantities at transparency. Finally the three different 
regimes of the afterglow are described within the fully radiative and radial approxima- 
tions: the ultrarelativistic, the relativistic and the nonrelativistic regimes. The best fit 
of the theory leads to an unequivocal identification of the "long GRBs" as extended 
emission occurring at the afterglow peak (E-APE). The relative intensities, the time 
separation and the hardness ratio of the P-GRB and the E-APE are used as distinctive 
observational test of the EMBH theory and the excellent agreement between our theoret- 
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ical predictions and the observations arc documented. The afterglow power-law indexes 
in the EMBH theory are compared and contrasted with the ones in the literature, and 
no beaming process is found for GRB 991216. Finally, some preliminary results relating 
the observed time variability of the E-APE to the inhomogeneities in the interstellar 
medium are presented, as well as some general considerations on the EMBH formation. 
The issue of the GSTS paradigm will be the object of a forthcoming publication and 
the relevance of the iron-lines observed in GRB 991216 is shortly reviewed. The gen- 
eral conclusions are then presented based on the three fundamental parameters of the 
EMBH theory: the dyadosphere energy, the baryonic mass of the remnant, the interstel- 
lar medium density. An in depth discussion and comparison of the EMBH theory with 
alternative theories is presented as well as indications of further developments beyond 
the radial approximation, which will be the subject of paper II in this series0 Future 
needs for specific GRB observations are outlined. 

Keywords: Afterglow, electromagnetic black hole theory, gamma-ray bursts 

1. Introduction 

1.1. The physical and astrophysical background 

Gamma-ray bursts (GRBs) are rapidly fuelling one of the broadest scientific pursuit 
in the entire field of science, both in the observational and theoretical domains. 
Following the discovery of GRBs by the Vela satellites,Elthe observations from the 
Compton satellite and BATSE^ had shown the isotropic distribution of the GRBs 
strongly suggesting a cosmological nature for their origin. It was still through the 
data of BATSE that the existence of two families of bursts, the "short bursts" and 
the "long bursts" was presented, opening an intense scientific dialogue on their 



An enormous momentum was gained in this field by the discovery of the after- 
glow phenomena by the BeppoSAX satellite and the optical identification of GRBs 
which have allowed the unequivocal identification of their sources at cosmological 
distances.^ It has become apparent that fluxes of 10^'* erg/s are reached: during 
the peak emission the energy of a single GRB equals the energy emitted by all the 
stars of the Universe.l^ 

From an observational point of view, an unprecedented campaign of observations 
is at work using the largest deployment of observational techniques from space with 
the sateUites CGRO-BATSE, Beppo-SAX*^, Chandra^ R-XTE'i, XMM-Newton°, 
HETE-2f, as well as the HST^, and from the ground with optical (KECK'^, VLT') 
and radio (VLAj) observatories. The further possibility of examining correlations 



^See http://space.mit.edu/HETE/" 
^See http://www.stsci.edu/ 
'^See http:/ /www2.keck.hawaii.edu:3636/ 
'See http:/ /www. eso.org/projects/vlt/ 
JSee http:/ /www. aoc.nrao.edu/vla/html/VLAhomc.shtml 
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with the detection of ultra high energy cosmic rays, UHECR for short, and in 
coincidence neutrinos should be reachable in the near future thanks to developments 
of AUGER'^ and AMANDA' (see also Halzen, 200(}2>). 



E-APE 




Fig. 1. Selected events in the EMBH theory are represented. For each one the values of the 
energy density of the medium and the distances from the EMBH, in the laboratory frame and in 
logarithmic scale, are given. 

From a theoretical point of view, GRBs offer comparable opportunities to de- 
velop entire new domains in yet untested directions of fundamental science. For the 
first time within the theory based on the vacuum polarization process occurring in 
an electromagnetic black hole, the EMBH theory, see Fig.^ the opportunity exists 
to theoretically approach the following fundamental issues: 

(1) The extremely relativistic hydrodynamic phenomena of an electron-positron 
plasma expanding with sharply varying gamma factors in the range 10^ to 10'* 
and the analysis of the very high energy collision of such an expanding plasma 
with baryonic matter reaching intensities 10'^^ larger than the ones usually 
obtained in Earth-based accelerators. 

(2) The bulk process of vacuum polarisati on cr eated by over criti cal electromagnetic 
fields, in the sense of Heisenberg, EuleJ^and SchwingeJ^. This longly sought 
quantum ultrarclativistic effect has not beei i yet unequivocally observed in 
heavy ion collision on the Earth! I ' ■ ^ 11 I The difficulty of the heavy ion 
collision experiments appears to be that the overcritical field is reached only for 

''See http://www.auger.org/ 

'See http:/ /amanda. berkeley.edu/amanda/amanda. html 
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time scales of the order h/mpC^, which is much shorter than the characteristic 
time for the e^e^ pair creation process which is of the order of h/rrieC^, where 
nip and me are respectively the proton and the electron mass. It is therefore 
very possible that the first appearance of such an effect occurs in the strong 
electromagnetic fields developed in astrophysical conditions during the process 
of gravitational collapse to an EMBH, where no problem of confinement exists. 
(3) A novel form of energy source: the extractable energy of a black hole. The 
enormous energies released almost instantly in the observed GRBs, points to the 
possibility that for the first time we are witnessing the release of the extractable 
energy of an EMBH, during the process of gravitational collapse itself. We can 
compute and have the opportunity to study all general relativistic as well as the 
associated ultrahigh energy quantum phenomena as the horizon of the EMBH 
is approached and is being formed. 

It is clear that in approaching such a vast new field of research, implying pre- 
viously unobserved relativistic regimes, it is not possible to proceed as usual with 
an uncritical comparison of observational data to theoretical models within the 
classical schemes of astronomy and astrophysics. Some insight to the new approach 
needed can be gained from past experience in the interpretation of relativistic effects 
in high energy particle physics as well as from the explanation of some observed 
relativistic effects in the astrophysical domain. Those relativistic regimes, both in 
physics and astrophysics, are however much less extreme than those encountered 
in GRBs. 

There are three major new features in relativistic systems which have to be 
properly taken into account: 

(1) Practically all data on astronomical and astr ophysical systems is acquired by 
using photon arrival times. It was EinsteirU^ at the very initial steps of special 
relativity who cautioned about the use of such an arrival time analysis and 
stated that when dealing with objects in motion proper care should be taken in 
defining the time synchronization procedure in order to construct the correct 
space-time coordinate grid (see Fig.|21). It is not surprising that as soon as the 
first relativistic bulk motion effects were observed their interpretations within 
the classical framework of astrophysics led to the concept of "superluminal" 
motio n. T hese were observations of extragalactic radio sources, with gamma 
fact or 10 and of microquasars in our own galaxy with gamma 

factonCl 

~ 5. It has been recognizecM^ that no "superluminal" motion exists if the 
prescriptions indicated by Einstein are used in order to establish the correct 
space-time grid for the astrophysical systems. In the present context of GRBs, 
where the gamma factor can easily surpass 10^, the direct application of clas- 
sical concepts leads to enormous "superluminal" behaviours (see Tab. An 
approach based on classical arrival time considerations as sometimes done in 
the current literature completely subverts the causal relation in the observed 



On the structure of the burst and afterglow of Gamma-Ray Bursts I: the radial approximation 5 



R 







R 



Fig. 2. This qualitative diagram illustrates the relation between the laboratory time interval At 
and the arrival time interval Ata for a pulse moving with velocity v in the laboratory time (solid 
line). We have indicated here the case where the motion of the source has a nonzero acceleration. 

The arrival time is measured using light signals emitted by the pulse (dotted lines). Ro is the 
distance of the observer from the EMBH, to is the laboratory time corresponding to the onset of 
the gravitational collapse, and r is the radius of the expanding pulse at a time t = to + At.^ 



astrophysical phenomenon. 
(2) One of the clear successes of rclativistic field theories has been the understand- 
ing of the role of four-momentum conservation laws in multiparticle collisions 
and decays such as in the reaction: n — > p+e~+Pe- From the works of Pauli and 
Fermi it became clear how in such a process, contrary to the case of classical 
mechanics, it is impossible to analyze a single term of the decay, the electron 
or the proton or the neutrino or the neutron, out of the context of the global 
point of view of the relativistic conservation of the total four momentum of 
the system. This in turn involves the knowledge of the system during the entire 
decay process. These rules are routinely used by workers in high energy particle 
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physics and have become part of their cultural background. If we apply these 
same rules to the case of the relativistic system of a GRB it is clear that it is just 
impossible to consider a part of the system, e.g. the afterglow, without taking 
into account the general conservation laws and whole relativistic history of the 
entire system. The description of the afterglow alone, as has been given at times 
in the literature, indeed possible within the framework of classical astronomy 
and astrophysics, is not viable in a relativistic astrophysics context where the 
space-time grid necessary for the description of the afterglow depends on the 
entire previous relativistic part of the worldline of the system (see also section 
HH). 

(3) The lifetime of a process has not an absolute meaning as special and general 
relativity have shown. It depends both on the inertial reference frame of the lab- 
oratory and of the observer and on their relative motion. Such a phenomenon, 
generally expressed in the "twin paradox" , has been extensively checked and 
confirmed to extremely h igh accuracy as a byproduct of the elementary particle 
physics (g-2) experiment.!^ This situation is much more extreme in GRBs due 
to the very large (in the range 10^-10'*) and time varying (on time scales rang- 
ing from fractions of seconds to months) gamma factors between the comoving 
frame and the far away observer (see Fig. Moreover in the GRB context 
such an observer is also affected by the cosmological recession velocities of its 
local Lorentz frame. 

1.2. The Relative Space-Time Transformations: the RSTT 
paradigm and current scientific literature 

Here are some of the reasons why we have recently presented a basic relative space- 
time transformation (RSTT) paradig nPl to be applied prior to the interpretation 
of GRB data. 

The first step is the establishment of the governing equations relating: 

a) The comoving time of the pulse (r) 

b) The laboratory time (t) 

c) The arrival time at the detector (ta) 

d) The arrival time at the detector corrected for cosmological expansion {f^) 

The book-keeping of the four different times and corresponding space variables must 
be done carefully in order to keep the correct causal relation in the time sequence 
of the events involved. 

As formulated the RSTT paradigm contains two parts: the first one is a nec- 
essary condition, the second one a sufficient condition. The first part reads: "the 
necessary condition in order to interpret the GRB data, given in terms of the arrival 
time at the detector, is the knowledge of the entire worldline of the source from the 
gravitational collapse" . 

Clearly such an approach is in contrast with articles in the current literature 
which emphasize either some qualitative description of the sources or some quanti- 
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tative description of the afterglow era by itself. 

In the current literature several attempts have addressed the is sue of the so urces 
of GRBs. They include scenarios of binary neutron stars mergers, "^3E32SI^ black 
hole / white dwari^'' and black hole / neutron star b inaries, ^21 hypernovaepSI 
failed supernovae or coUapsars ; 29^313 g^pj-^novae!^^ Only those based on binary 
neutron stars have reached the stage of a definite model and detailed quantitative 
estimates have been made. In this case, however, various problems have surfaced: in 
the general energetic s whic h cannot be greater than 3 x 10^^ erg, in the explana- 
tion of "long bursts" j'^^^ l ^^^ l and in the observed location of the GRB sources in star 
forming regionsP^ In the remaining cases attention was directed to a qualitative 
analysis of the sources without addressing the overall problem from the source to 
the observations. The necessary details to formulate the equations of the dynamical 
evolution of the system are generally missing. 

Other models in the literature have addressed the problem of only fitting the 
data of the afterglow observations by a phenomenological analysis. They are sepa- 
rated into two major classes: 

The "internal shock model", first introduced by Rees & Meszaros (1994), ESI 
by far the most popula r on e, has been developed in many differe nt as pects, e.g. by 
Paczynski & Xu (1994)E3sari & Piran (1997),^^ Fenimore (1999^31 and Fenimore 
et al. (1999j3ni. The underlying assumption is that all the variabilities of GRBs in 
the range Ai ~ 1 ms up to the overall duration T of the order of 50 s are determined 
by a yet undetermined "inner engine" . The difficulties of explaining the long time 
scale bursts by a single explosive model has evolved into a subclass of approaches 
assuming an "inner engine" with extended activity (see e.g. Piran, 2001,1^ and 
references therein). 

The "external shock model", also introduced by Meszaros & Rees (1993),E21 
is less popular today. It relates the GRB light curves and time variabilities to 
interactions of a single thin blast wave with clouds in the external medium. The 
interesting possibility has been recognized within this model, that GRB light curves 
"are tomographic images of the density distribution of the medium surrounding the 
sources of GRBs" (Dermer & Mitman, 199^321) 

see also Dermer, Chiang & Bottcher 

( 1999) Dermer (2002^33 

and references therein. In this case, the structure of the 
burst is assumed not to depend directly on the "inner engine" (see e.g. Piran, 
2001,1^ and references therein). 

All these works encounter the above mentioned difficulty: they present either 
a purely qualitative or phenomenological or a piecewise description of the GRB 
phenomenon. By neglecting the earlier phases, their space-time grid is undefined 
and as we will explicitly show in the following, results are reached at variance from 
the ones obtained in a complete and unified description of the GRB phenomenon. 
We show in the following how such a unified description naturally leads to new 
characteristic features both in the burst and afterglow of GRBs. 
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1.3. The EMBH Theory 

In a series of papers, we have developed the EMBH theorjS^ which has the ad- 
vantage, despite its simplicity, that all eras following the process of gravitational 
collapse are described by precise field equations which can then be numerically 
integrated. 

Starting from the vacuum polarization process a la Heisenberg-Euler-SchwingeJ^^^ 
in the overcritical field of an EMBH first comi 3uted in Damour & Ruffini (1975),E3 
we have developed the dyadosphere concept l^sl 

The dynamics of the e^e~-pairs and electromagnetic radiation of the plasma 
generated in the dyadosphere propagating away from the EMBH in a sharp pulse 
(PEM pulse) has been studied by the Rome group and validated by the numerical 
codes developed at Livermore Lab.*^^ 

The collision of the still optically thick e+e~-pairs and electromagnetic radiation 
plasma with the baryonic matter of the remnant of the progenitor star h as b een 
again studied by the Rome group and validated by the Livermore Lab codesl^The 
further evolution of the sharp pulse of pairs, electromagnetic radiation and baryons 
(PEMB pulse) has been followed for increasing values of the gamma factor until 
the condition of transparency is reached!^ 

As this PEMB pulse reaches transparency the proper GRB (P-GRB) is emitteci^l 
and a pulse of accelerated baryonic matter (the ABM pulse) is injected into the 
interstellar medium (ISM) giving rise to the afterglow. 



1.4. The GRB 991216 as a prototypical source 

Until this stage, the EMBH theory has been done from first principles based on the 
exact solutions of the Einstein-Maxwell equations implied by the EMBH uniqueness 
theorem as well as on the quantum description of the vacuum polarization process 
in overcritical electromagnetic fields. Turning now to the afterglow, the variety of 
physical situations that can possibly be encountered are very large and far from 
unique: the description from first principles is just impossible. We have therefore 
proceeded to properly identify what we consider a prototypical GRB source and 
to develop a theoretical framework in close correspondence with the observational 
data. 

We present the criteria which have guided us in the selection of the GRB source 
to be used as a prototype before proceeding to an uncritical comparison with 
the theory. It is now clear, since the observations of GRB 980425, GRB 991216, 
GRB 970514 and GRB 980326 that the afterglow phenomena can present, espe- 
cially in the optical and radio wavelengths, features originating from phenomena 
spatially and causally distinct from the GRB phenomena. There is the distinct 
possibility that phenomena related to a supernova can be erroneously attributed 
to a GRB. This problem has been clearly addressed by the GRB supernova time 
sequence (GSTS) paradigm in which the time sequence of the events in the GRB 
supernova phenomena has been outlined.'^ This has led to the novel concept of an 
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induced supernova.'-' This problem will be addressed in a forthcoming paper 
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Fig. 3. a) The peak emission of GRB 991216 as seen by BATSE (Reproduced from BATSE 
Rapid Burst Response^^); b) The afterglow emission of GRB 991216 as seen by XTE and Chandra 
(reproduced from Halpern et al., 2000^'*) 



In view of these considerations we have selected GRB 991216 as a prototypical 
case (see Fig.Ol for the following reasons: 
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(1) GRB 991216 is one of the strongest GRBs in X-rays and is also quite general in 
the sense that it shows relevant cosmological effects. It radiates mainly in X-ra ys 
and in 7-rays and less than 3% is emitted in the optical and radio bands!^ 

(2) The excellent data obtained by BATSE on the bursl^ is complemented by the 
data on the afterglow acquired by ChandrJ^ and RXTeEBIaIso superb data 
have been obtained from spectroscopy of the iron lines!^ 

(3) A value for the slope of the energy er niss ion during the afterglow as a function 
of time has been obtained: n = —1.64^ and n = — 1.616 ± 0.067.*^ 



1.5. The interpretation of the burst structure: the IBS paradigm 
and the different eras of the EMBH theory 

The comparison of the EMBH theory with the data of the GRB 991216 and its 
afterglow has naturally led to a new paradigm for the interpretation of the burst 
structures (IBS paradigm)) of GRBs. The IBS paradigm reads: "In GRBs we 
can distinguish an injector phase and a beam-target phase. The injector phase in- 
cludes the process of gravitational collapse, the formation of the dyadosphere, as 
well as Era I (the PEM pulse). Era II (the engulfment of the baryonic matter of the 
remnant) and Era III (the PEMB pulse). The injector phase terminates with the P- 
GRB emission. The beam-target phase addresses the interaction of the ABM pulse, 
namely the beam generated during the injection phase, with the ISM as the target. 
It gives rise to the E-APE and the decaying part of the afterglow". The detailed 
presentations of these results are the main topic of this article. 

We recall that the injector phase starts from the moment of gravitational 
collapse and encompasses the following eras: 

The Zeroth Era: the formation of the dyadosphere. In section [21 we review the 
basic scientific results which lie at the basis of the EMBH theory: the black hole 
uniqueness theorem, the mass formula of an EMBH, the process of vacuum polar- 
ization in the field of an EMBH. We also point out how after the discovery of the 
GRB afterglow the reexamination of these results has led to the novel concept of 
the dyadosphere of an EMBH. We have investigated this concept in the simplest 
possible case of an EMBH depending only on two parameters: the mass and charge, 
corresponding to the Reissner-Nordstrom spacetime. We recall the definition of the 
energy Edya of the dyadosphere as well as the spatial distribution and energetics of 
the e+e~ pairs. See Fig.^ 

In order to analyse the time evolution of the dyadosphere we give in the three 
following sections the theoretical background for the needed equations. 

In section |31 we give the general relativistic equations governing the hydrody- 
namics and the rate equations for the plasma of e'*'e^-pairs. 

In section0]we give the governing equations relating the comoving time r to the 
laboratory time t corresponding to an incrtial reference frame in which the EMBH 
is at rest and finally to the time measured at the detector ta which, to finally get 
if, must be corrected to take into account the cosmological expansion. 



On the structure of the burst and afterglow of Gamma- Ray Bursts I: the radial approximation 11 




plasma (e*er y) 



Fig. 4. The dyadosphere of a Reissner-Nordstrom black hole can be represented as constituted 
by a concentric set of shells of capacitors, each one of thickness h/rrieC and producing a number 
of e+e" pairs of the order of ~ Q/e on a time scale of 10~^^ s, where Q is the EMBH charge. 
The shells extend in a region Ar, from the horizon r+ to the dyadosphere outer radius r^s (see 
text). The system evolves to a thermalised plasma configuration. 

In sectional we describe the numerical integration of the hydrodynamical equa- 
tions and the rate equation developed by the Rome and Livermore groups. This 
entire research program could never have materialized without the fortunate inter- 
action between the complementary computational techniques developed by these 
two groups. The validation of the results of the Rome group by the fully general 
relativistic Livermore codes has been essential both from the point of view of the 
validity of the numerical results and the interpretation of the scientific content of 
the results. 

The Era I: the PEM pulse. In section |31 by the direct comparison of the inte- 
grations performed with the Rome and Livermore codes we show that among all 
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possible geometries the e+e" plasma moves outward from the EMBH reaching a 
very unique relativistic configuration: the plasma self-organizes in a sharp pulse 
which expands in the comoving frame exactly by the amount which compensates 
for the Lorentz contraction in the laboratory frame. The sharp pulse remains of 
constant thickness in the laboratory frame and self-propels outwards reaching ul- 
trarelativistic regimes, with gamma factors larger than 10^, in a few dyadosphere 
crossing times. We recall that, in analogy with the electromagnetic (EM) pulse 
observed in a thermonuclear explosion on the Earth, we have defined this more 
energetic pulse formed of electron-positron pairs and electromagnetic radiation a 
pair-electromagnetic-pulse or PEM pulse. 

The Era II: We describe the interaction of the PEM pulse with the baryonic 
remnant of mass Mb left over from the gravitational collapse of the progenitor 
star. We give the details of the decrease of the gamma factor and the corresponding 
increase in the internal energy during the collision. The dimensionless parameter 
B = Mbc^ / Edya which measures the baryonic mass of the remnant in units of the 
Edya is introduced. This is the second fundamental free parameter of the EMBH 
theory. 

The Era III: We describe in section|Slthe further expansion of the e+e^ plasma, 
after the engulfment of the baryonic remnant of the progenitor star. By direct 
comparison of the results of integration obtained with the Rome and the Livermore 
codes it is shown how the pair-electromagnetic-baryon (PEMB) plasma further 
expands and self organizes in a sharp pulse of constant length in the laboratory 
frame (see Fig. O. We have examined the formation of this PEMB pulse in a wide 
range of values 10^^ < B < 10^^ of the parameter B, the upper limit corresponding 
to the limit of validity of the theoretical framework developed. 

In section it is shown how the effect of baryonic matter of the remnant, 
expressed by the parameter B, is to smear out all the detailed information on 
the EMBH parameters. The evolution of the PEMB pulse is shown to depend only 
on Edya and B: the PEMB pulse is degenerate in the mass and charge parameters 
of the EMBH and rather independent of the exact location of the baryonic matter 
of the remnant. 

In section El the relevant thermodynamical quantities of the PEMB pulse, the 
temperature in the different frames and the e^e~ pair densities, are given and the 
approach to the transparency condition is examined. Particular attention is given 
to the gradual transfer of the energy of the dyadosphere Edya to the kinetic energy 
of the baryons Esaryons during the optically thick part of the PEMB pulse. 

In section 1111 as the condition of transparency is reached, the injector phase 
is concluded with the emission of a sharp burst of electromagnetic radiation and 
an accelerated beam of highly relativistic baryons. We recall that we have respec- 
tively defined the radiation burst (the proper GRB or for short P-GRB) and the 
accelerated-baryonic-matter (ABM) pulse. By computing for a fixed value of the 
EMBH different PEMB pulses corresponding to selected values of B in the range 
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Fig. 5. Comparison of gamma factor for the one-dimensional (1-D) hydrodynamic calculations 
(Livermore code) and slab calculations (Rome code) as a function of the radial coordinate (in units 
of dyadosphere radius) in the laboratory frame. The calculations show an excellent agreement. 



[lO~^-10^^], it has been possible to obtain a crucial universal diagram which is 
reproduced in Fig^l In the limit of B ^ 10~^ or smaller almost all E^ya is emitted 
in the P-GRB and a negligible fraction is emitted in the kinetic energy EBaryons 
of the baryonic matter and therefore in the afterglow. On the other hand in the 
limit B 10~^ which is also the limit of validity of our theoretical framework, 
almost all Edya is transferred to Esaryons and gives origin to the afterglow and the 
intensity of the P-GRB correspondingly decreases. We have identified the limiting 
case of negligible values of B with the process of emission of the so called "short 
bursts" . A complementary result reinforcing such an identification comes from the 
thermodynamical properties of the P-GRB: the hardness of the spectrum decreases 
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Fig. 6. At the transparent point, the energy radiated in the P-GRB (the solid line) and the final 
kinetic energy Esaryons of baryonic matter (the dashed line) in units of the total energy of the 
dyadosphere (Sdya) a^re plotted as functions of the B parameter. 

for increasing values of B, see Fig. 13 

The injector phase is concluded by the emission of the P-GRB and the ABM 
pulse, as the condition of transparency is reached. 

The beam-target phase, in which the accelerated baryonic matter (ABM) 
generated in the injector phase collides with the ISM, gives origin to the afterglow. 
Again for simplicity we have adopted a minimum set of assumptions: 

(1) The ABM pulse is assumed to collide with a constant homogeneous interstellar 
medium of number density fiism ~ Icm^"^. The energy emitted in the collision 
is assumed to be instantaneously radiated away (fully radiative condition) . The 
description of the collision and emission process is done using spherical sym- 
metry, taking only the radial approximation neglecting all the delayed emission 
due to off-axis scattered radiation. 

(2) Special attention is given to numerically compute the power of the afterglow 
as a function of the arrival time using the correct governing equations for the 
space-time transformations in line with the RSTT paradigm. 

(3) Finally some approximate solutions are adopted in order to determine the power 
law exponents of the afterglow flux and compare and contrast them with the 
observational results as well as with the alternative results in the literature. 
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Fig. 7. The energy corresponding to the peak of the photon number spectrum in the P-GRB as 
measured in the laboratory frame is plotted as function of the B parameter. 



In this paper we only consider the above mentioned radial approximation and a 
spherically symmetric distribution in order to concentrate on the role of the correct 
space-time transformations in the RSTT paradigm and illustrate their impact on 
the determination of the power law index of the afterglow. This topic has been 
seriously neglected in the literature. Details of the role of bea ming and on the 
diffusion due to off-axis emission will be studied elsewhere .l^^^^ 

We can now turn to the two eras of the beam-target phase: 

The Era IV: the ultrarelativistic and relativistic regimes in the afterglow. In 
section 1121 the hydrodynamic relativistic equations governing the collision of the 
ABM pulse with the interstellar matter are given in the form of a set of finite dif- 
ference equations to be numerically integrated. Expressions for the internal energy 
developed in the collision as well as for the gamma factor are given as a function 
of the mass of the swept up interstellar material and of the initial conditions. In 
section 1171 the infinitesimal limit of these equations is given as well as analytic 
power-law expansions in selected regimes. 

The Era V: the approach to the nonrelativistic regimes in the afterglow. In 
section [T51 it is stressed that this last era often discussed in the current literature 
can be described by the same equations used for era IV. 

Having established all the governing equations for all the eras of the EMBH 
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theory, we can proceed to compare and contrast the predictions of this theory with 
the observational data. 

1.6. The Best fit of the EMBH theory to the GRB 991216: the 
global features of the solution 

As expressed in section we have proceeded to the identification of the only two 
free parameters of the EMBH theory, Edya and B, by fitting the observational data 
from R-XTE and Chandra on the decaying part of the GRB 991216 afterglow. 
The afterglow appears to have three different parts: in the first part the luminosity 
increases as a function of the arrival time, it then reaches a maximum and finally 
monotonically decreases. In Fig. |H1 we show how such a fit is actually made and 
how changing the two free parameters affects the intensity and the location in time 
of the peak of the afterglow. The best fit is obtained for Edtja — 4.83 x 10^^ erg and 
S = 3 X 10-3. 

Having determined the two free parameters of the theory, we have integrated the 
governing equations corresponding to these values and then obtained for the first 
time the complete history of the gamma factor from the moment of gravitational 
collapse to the latest phases of the afterglow observations (see Fig.Ell. We have also 
determined the different regimes encountered in the relation between the laboratory 
time and the detector arrival time within the RSTT paradigm (see Fig. Il()|l . We 
have thus determined the entire space-time grid of the GRB 991216 by giving 
(see Tab. ^ the radial coordinate of the GRB phenomenon as a function of the 
four coordinate time variables. A quick glance to Tab. Q shows how the extreme 
relativistic regimes at work lead to enormous superluminal behaviour (up to lO^c!) 
if the classical astrophysical concepts are adopted using the arrival time as the 
independent variable. In turn this implies that any causal relation based on classical 
astrophysics and the arrival time data, as often found in the current GRB literature, 
is incorrect. 

1.7. The explanation of the "long hursts" and the identification of 
the proper gamma ray burst(P-GRB) 

In sectional having determined the two free parameters of the EMBH theory, we 
analyze the theoretical predictions of this theory for the general structure of GRBs. 
The first striking result, illustrated in Fig.^] shows that the peak of the afterglow 
emission coincides both in intensity and in arrival time (19.87 s) with the average 
emission of the long burst observed by BATSE. For this we have introduced the 
new concept of extended afterglow peak emission (E-APE). Once the proper space- 
time grid is given (see Tab. ^ it is immediately clear that the E-APE is generated 
at distances of 5 x 10^^ cm from the EMBH. The long bursts are then identified 
with the E-APEs and are not bursts at all: they have been interpreted as bursts 
only because of the high threshold of the BATSE detectors (see Fig. Ill|l . Thus the 
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Fig. 8. a) Afterglow luminosity computed for an EMBH of E^ya = 5.29 X 10®^ erg, E^i/a = 
4.83 X 10^3 erg, E^y^ = 4.49 X 10^5 erg and B = 3 X IQ-^. b) for the E^y^ = 4.83 X 10^3, -we 
give the afterglow luminosities corresponding respectively to B = 9 X 10"'^, 6 X 10"'^, 3 X 10"^, 
1 X 10-3, 7 X 10"*, 4 X 10"*. 



long standing unsolved problem of explaining the long is radically 

resolved. 

Still in sectional the search for the identification of the P-GRB in the BATSE 
data is described. This identification is made using the two fundamental diagrams 
shown in Fig. 1121 and Fig.^] Having established the value of Edya — 4.83 x 10^"^ erg 
and of i? = 3 X lO"'^, it is possible from the dashed line and the solid line in Fig. 1121 
to evaluate the ratio of the energy Ep-crb emitted in the P-GRB to the energy 
Esaryons emitted in the afterglow corresponding to the determined value of B, see 
the vertical line in Fig. ^1 We obtain Ep-crb / Esaryons = 1-58 x 10"^, which 
gives Ep-GRB = 7.54 x 10*^^ erg. Having so determined the theoretically expected 
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Fig. 9. The theoretically computed gamma factor for the parameter values -B^i/a = 4.83 X 10 
erg, S = 3 X 10~^ is given as a function of the radial coordinate in the laboratory frame. The 
corresponding values in the comoving time, laboratory time and arrival time are given in Tab. 
The different eras indicated by roman numerals are illustrated in the text (see sections l6l7l8ll2ll3l , 
while the points 1,2,3,4,5 mark the beginning and end of each of these eras. The points Pi^ and 
Pa mark the maximum of the afterglow flux, respectively in emission time and in arrival time'^ 
(see and sections ll2ll7i . The point 6 is the beginning of Phase D in era V (see sections ll3ll7i . At 
point 4 the transparency condition is reached and the P-GRB is emitted. 

intensity of the P-GRB, a second fundamental observable parameter, which is also 
a function of Edya and B, is the arrival time delay between the P-GRB and the 
peak E-APE, determined in Fig.^1 From Tab.^ we have that the detector arrival 
time of the P-GRB occurs at 8.41 x 10^^ s, corresponding to a radial coordinate of 
1.94 X lO'^* cm, a comoving time of 21.57 s, a laboratory time of 6.48 x 10"^ s and an 
arrival time of 4.21 x 10^^ s. At this point, the gamma factor is 310.1. The peak of 
the E-APE occurs at a detector arrival time of 19.87 s, corresponding to a radial 
coordinate of 5.18 x IQ^^ cm, a comoving time of 5.85 x 10^ s, a laboratory time 
of 1.73 X 10^ s and an arrival time of 9.93 s (see Tab. QJ. The delay between the 
P-GRB and the peak of the E-APE is therefore 19.78 s, see Fig. [HI The theoretical 
prediction on the intensity and the arrival time uniquely identifies the P-GRB with 
the "precursor" in the GRB 991216 (see Fig. 13). Moreover, the hardness of the 
P-GRB spectra is also evaluated in this section. As pointed out in the conclusions, 
the fact that both the absolute and relative intensities of the P-GRB and E-APE 
have been predicted within a few percent accuracy as well as the fact that their 
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Fig. 10. Relation between the arrival time (if) measured at the detector and the laboratory time 
(t) measured at the GRB source. The solid curve is computed using the exact formula given in 
Eq.J^. The dashed-dotted curve is computed using the approximate formula given in EQ. I39i 
and often used in the current literature. We distinguish four different phases. Phase A: There is 
a linear relation between t and tjj, given by EQ. I134t in the text (dashed line). Phase B; There 
is an "effective" power-law relation between t and if, given by Eg. 1139^ (dotted line). Phase 
C: No analytic formula holds and the relation between t and tjj has to be directly computed by 
the integration of the complete equations of energy and momentum conservation fEas. llU4llU5l ). 
Phase D; As the gamma factor approaches 7=1, the relation between t and 4jJ asymptotically 
goes to t = (light gray line).-'^ 

arrival time has been computed with the precision of a few tenths of miUiseconds, 
see Tab. Hand Fig. ^1 can be considered one of the major successes of the EMBH 
theory. 



1.8. On the power-laws, beaming and temporal structures in the 
afterglow of GRB 991216. 

In section lTTI a piecewise description of the afterglow b y th e expansion of the funda- 
mental hydrodynamical equations given by Taub (1948 and Landau & LifshitJHIl 
have allowed the determination of a power-law index for the dependence of the af- 
terglow luminosity on the photon arrival time at the detector. It is evident that 
the determination of the power-law index is very sensitive to the basic assumptions 
made for the description of the afterglow, as well as to the relations between the 
different temporal coordinates which have been clarified by the RSTT paradigmP 
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Table 1. Gamma factors for selected events and their space-time coordinates. The points marked 
l,2,3,4,5,6,Pi,Pyi are the same reported in Fig. |5] while the point F is the endpoint of the sim- 
ulation. It is particularly important to read the last column, where the apparent motion in the 
radial coordinate, evaluated in the arrival time at the detector, leads to an enormous "superlumi- 
nal" behaviour, up to 9.55 X 10* c. This illustrates well the impossibility of using such a classical 
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The different power-faw indices obtained are compared and contrasted witii tiie 
ones in tlie current iiterature (see Tab. 13) • As a byproduct of tfiis analysis, see also 
the conclusions, there is a perfect agreement between the observational data and 
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Fig. 11. Best fit of the afterglow data of Chandra, RXTE as well as of the range of variability of 
the BATSE data on the major burst, by a unique afterglow curve leading to the parameter values 
Edya = 4.83 X Vf'^erg,B = 3 X 10~^. The horizontal dotted fine indicates the BATSE noise 
threshold. On the left axis the luminosity is given in units of the energy emitted at the source, 
while the right axis gives the flux as received by the detectors. 



the theoretical predictions, implying that the assumptions adopted for the descrip- 
tion of the afterglow are valid and therefore that there is no evidence for a beamed 
emission in GRB 991216. 

In section El the role of the inhomogeneities in the interstellar matter has been 
analysed in order to explain the observed temporal substructures in the BATSE 
data on GRB 991216. From the data of Tab. ^ and the highly "superluminal" 
behaviour of the source in the region of the E-APE, it is concluded that the observed 
time variability in the intensity of the emission (A///) ~ 5 can be traced to 
inhomogeneities in the interstellar matter: (Ariism / nism) ~ 5- The typical size of 
the scattering region is estimated to be 5 x 10^^ cm, and these are the typical 
sizes and density contrasts found in interstellar clouds. Since the emission of the 
E-APE occurs at typical dimensions of the order of 5 x 10^^ cm, the observed 
inhomogeneities are probing the structure of the interstellar medium, and have 
nothing to do with the "inner engine" of the source. These conclusions, reached in 
the radial approximation of the afterglow adopted in this article, have been proved 
to hold in the more general case when off-radial emission is taken into account .E^^ 
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Fig. 12. Relative intensities of the E-APE (dashed line) and the P-GRB (solid line), as predicted 
by the EMBH theory corresponding to the values of the parameters determined in Fig. 1111 as 
a function of B, Details arc given in section [THl The vertical line corresponds to the value B = 
3 X 10"^. 

1.9. The observation of the iron lines in GRB 991216: on a 
possible GRB-supernova time sequence 

In section 1201 the program of using GRBs to further explore the region surround- 
ing the newly formed E MBH is carried one step further by using the observations 
of the emitted iron lines.l^ This gives us the opportunity to introduce the GRB- 
supernova time sequence (GSTS) paradigm and to introduce as well the novel con- 
cept of an induced supernova explosion. The GSTS paradigm reads: A massive 
GRB-progenitor star Pi of mass Mi undergoes gravitational collapse to an EMBH. 
During this process a dyadosphere is formed and subsequently the P-GRB and the 
E-APE are generated in sequence. They propagate and impact, with their photon 
and neutrino components, on a second supernova-progenitor star P2 of mass M2. 
Assuming that both stars were generated approximately at the same time, we expect 
to have AI2 < Mi. Under some special conditions of the thermonuclear evolution of 
the supernova-progenitor star P2, the collision of the P-GRB and the E-APE with 
the star P2 can induce its supernova explosion. 

Using the result presented in Tab. ^a-nd in a-U preceding sections, the GSTS 
paradigm is illustrated in the case of GRB 991216. Some general considerations on 
the nature of the supernova progenitor star are also advanced. 
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Fig. 13. The arrival time delay between the P-GRB and the peak of the E-APE is plotted as a 
function of the B parameter for three selected values of E^iya- 

Some general considerations on the EMBH formation are presented in sectionl^ 
The general conclusions are presented in section 1^ 

The understanding of all these points has led to the formulation of the second 
part, namely the sufficient condition of the RSTT paradigm which reads: "the 
necessary condition in order to interpret the GRB data, given in terms of the arrival 
time at the detector, is the knowledge of the entire worldline of the source from the 
gravitational collapse. In order to meet this condition, given a proper theoretical 
description and the correct governing equations, it is sufficient to know the energy 
of the dyadosphere and the mass of the remnant of the progenitor star" . 

2. The zeroth era: the process of gravitational collapse and the 
formation of the dyadosphere 

We first recall the three theoretical results which lie at the basis of the EMBH 
theory. 

In 1971 in the article "Introducing the Black Hole ",Elthe theorem was advanced 
that the most general black hole is characterized uniquely by three independent 
parameters: the mass-energy M, the angular momentum L and the charge Q making 
it an EMBH. Such an ansatz, which came to be known as the "uniqueness theorem" 
has turned out to be one of the most difficult theorems to be proven in all of physics 
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Fig. 14. A qualitative diagram showing the full picture of the model, with both P-GRB and 
E-APE. 



and mathema tics . The progress in the proof has been authoritatively summarized by 
Carter (1997).ESlThe situation can be considered satisfactory from the point of view 
of the physical and astrophysical considerations. Nevertheless some fundamental 
mathematical and physical issues concerning the most general perturbation analysis 
of an EMBH are still the topic of active scientific discussion.l^ 

In 1971 it was shown that the energy extractable from an EMBH is governed 
by the mass-energy formulapS 



2p+ 



with 



where 



4LV) < 1, 



S = 47rp^ = ATT{r\ + 



) 



16^ ( ^ ] Af2 



(1) 



(2) 



(3) 



is the horizon surface area, Mj^ is the irreducible mass, r+ is the horizon radius 
and /9+ is the quasi-spheroidal cylindrical coordinate of the horizon evaluated at 
the equatorial plane. Extreme EMBHs satisfy the equality in Eq.Q. Up to 50% of 
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the mass-energy of an extreme EMBH can in princi ple be extracted by a special 
set of transformations: the reversible transformations!"^ 

In 1975, generalizing some previous results of Zaumen (1975^^ and Gibbons 
(1975)J^Damour & Ruffini (19 75'P3 

showed that the vacuum polarization process 

a la . 
than 



a la Heisenberg-Euler-Schwingei^^^^ created by an electric field of strength larger 



S. = ^ (4) 
he 

can indeed occur in the field of a Kerr-Newmann EMBH. Here TOe and e are respec- 
tively the mass and charge of the electron. There Damour and Ruffini considered an 
axially symmetric EMBH, due to the presence of rotation, and limited themselves 
to EMBH masses larger then the upper limit of a neutron star for astrophysical 
applications. They purposely avoided all complications of black holes with mass 
smaller then the dual electron mass of the electron ( m* — -pf^ ~ "'piancfc \ 

may lead to quantum evaporation processes P They pointed out that: 

(1) The vacuum polarization process can occur for an EMBH mass larger than the 
maximum critical mass for neutron stars all the way up to 7.2 x 10^ Mq. 

(2) The process of pair creation occurs on very short time scales, typically , and 
is an almost perfect reversible process, in the sense defined by Christodoulou- 
Rufhni, leading to a very efficient mechanism of extracting energy from an 
EMBH. 

(3) The energy generated by the energy extraction process of an EMBH was found 
to be of the order of lO'^'* erg, released almost instantaneously. They concluded 
at the time "this work naturally leads to a most simple model for the explanation 
of the recently discovered ^-ray bursts". 

After the discovery of the afterglow of GRBs and the determination of the cos- 
mological distance of their sources we noticed the coincidence between the theoret- 
ically predicted energetics and the observed ones in Damour & Rufhni (1975)4^ we 

returned to our theoretical results developing some new basic theoretical concepts j '^^ l ^^ l ^^ l ^^ l ^^ l 
which have led to the EMBH theory. 

As a first simplifying assumption we have developed our considerations in the 
absence of rotation with spherically symmetric distributions. The space-time is then 
described by the Reissner-Nordstrom geometry, whose spherically symmetric metric 
is given by 

d^s = gtt{r)d^t + grr{r)d^r + r'^d^O + sin^ ed^<t> , (5) 

where 5tt(r) = - 1 - 2GM + Oig ^ ^^'^{^r) and grr{r) ^ a-'^{r). 

The first new result we obtained is that the pair creation process does not occur 
at the horizon of the EMBH: it extends over the entire region outside the horizon in 
which the electric field exceeds the critical value given by Eq. 0| Since the electric 
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field in the Reissner-Nordstrom geometry has only a radial component given bjl^ 



this region extends from the horizon radius 

r+ = 1.47 • 10V(1 + -C^) cm 
out to an outer radiu^^ 



(6) 



(7) 



mc J \ c 
1.12 • 10^\//7C cm. 



) \m) \qpj V\/G, 



Q 

GM 



C8) 



M 



where we have introduced the dimensionless mass and charge parameters ^jl — 



Q 



< 1, see Fig.H 



The second new result has been to realize that the local number density of 
electron and positron pairs created in this region as a function of radius is given by 



-(r) 



Q 



K mc / 



1 



(9) 



47rr2 I rr:; ) e 

and consequently the total number of electron and positron pairs in this region is 

(r* 



N°+ 



1 



Tl 
mc 



(10) 



Scvl- 



where Qc 

The total number of pairs is larger by an enormous factor r* / [fi/mc) > 10^® 
than the value Q/e which a naive estimate of the discharge of the EMBH would 
have predicted. Due to this enormous amplification factor in the number of pairs 
created, the region between the horizon and r* is dominated by an essentially high 
density neutral plasma of electron-positron pairs. We have defined this region as the 
dyadosphere of the EMBH from the Greek duas, duadsos for pairs. Consequently we 
have called r* the dyadosphere radius r* = rds- ^^'*^^"^^' The vacuum polarization 
process occurs as if the entire dyadosphere are subdivided into a concentric set of 
shells of capacitors each of thickness h/nieC and each producing a number of e+e~ 
pairs on the order of ^ Q/e (see Fig.0J. The energy density of the electron-positron 
pairs is given by 



e{r) = 



1- ( — 



(11) 



(see Figs. 2-3 of Preparata, Ruffini & Xue, 1998cl^. The total energy of pairs 
converted from the static electric energy and deposited within the dyadosphere is 
then 

2' 



-fc'dya - T; ) 

2 r+ Tds 



I- I 

rds 



(12) 
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Fig. 15. The energy extracted by the process of vacuum polarization is plotted (solid lines) 
as a function of the mass M in solar mass units for selected values of the charge parameter 
5 = 1,0.1,0.01 (from top to bottom) for an EMBH, the case 5 = 1 reachable only as a limiting 
process. For comparison we have also plotted the maximum energy extractable from an EMBH 
(dotted lines) given by eq. 0. Details in Preparata, Ruffini & Xue (2001).^^ 



As we will see in the following this is one of the two fundamental parameters of 
the EMBH theory (see Fig.HB)). In the hmit ^ ^ 0, Eq.lO leads to E^y^ ^ 
which coincides with the energy extractable from EMBHs by reversible processes 
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(Mir — const.), namely Ebh — Mr = ^^p^see Fig. Due to the very large 
pair density given by Eq. @ and to the sizes of the cross-sections for the process 
e^e^ <^ 7 + 7, the system is expected to thcrmalize to a plasma configuration for 
which 



(13) 



where n°^^- is the total number density of e+e'-pairs created in the dyadosphere.l^^EHl 

The third new result which we have introduced for simplicity is that for a given 
Edya we have assumed either a constant average energy density over the entire 
dyadosphere volume, or a more compact configuration with energy density equal 
to the peak value. These are the two possible initial conditions for the evolution of 
the dyadosphere (see Fig. [T7jl . 
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Fig. 16. Selected lines corresponding to fixed values of the E^iya are given as a function of the 
two parameters /i 5, only the solutions below the continuous heavy line are physically relevant. The 
configurations above the continuous heavy lines correspond to unphysical solutions with r^s < ?"+ 



These three old and three new theoretical results permit a good estimate of the 
general energetics processes originating in the dyadosphere, assuming an already 
formed EMBH. In reality, if the data become accurate enough, the full dynamical 
description of the dyadosphere formation mentioned above will be needed in order 
to follow all the general relativistic effects and characteristic time scales of the 
approach to the EMBH horizoiJZ^»22iI4|75j^ g^jg^ section |5T] 
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Fig. 17. Two difFerent approximations for the energy density profile inside the dyadosphere. The 
first one (dashed fine) fixes the energy density equal to its peak value, and computes an "effective" 
dyadosphere radius accordingly. The second one (dotted line) fixes the dyadosphere radius to its 
correct value, and assumes an uniform energy density over the dyadosphere volume. The total 
energy in the dyadosphere is of course the same in both cases. The solid curve represents the real 
energy density profile. 



Below we shall concentrate on the dynamical evolution of the electron-positron 
plasma created in the dyadosphere. We shall first examine in the next three sections 
the governing equations necessary to approach such a dynamical description. 



3. The hydrodynamics and the rate equations for the plasma of 
e+e~-pairs 

The evolution of the e"^e~- pair plasma generated in the dyadosphere has been 
treated in two papersP^I^ We recall here the basic governing equations in the 
most general case in which the plasma fluid is composed of e~''e~ -pairs, photons 
and baryonic matter. The plasma is described by the stress-energy tensor 

T^"' ^pg^"' + {p + p)U^'U\ (14) 

where p and p are respectively the total proper energy density and pressure in the 
comoving frame of the plasma fluid and is its four- velocity, satisfying 



(15) 



30 R. Ruffini, C.L. Bianco, P. Chardonnet, F. Fraschetti, S.-S. Xue 



where and U*' are the radial and temporal contravariant components of the 
4- velocity. 

The conservation law for baryon number can be expressed in terms of the proper 
baryon number density ub 



- inBU'),t + ^ir^nBU'-)^r=0 . 



(16) 



The radial component of the energy-momentum conservation law of the plasma 
fluid reduces to 



dr dt 
-liP + P) 



dgr 



dr 



■([/■ 



r\2 







(17) 



The component of the energy-momentum conservation law of the plasma fluid equa- 
tion along a flow line is 

1 



P 







(18) 



Defining the total proper internal energy density e and the baryonic mass density 
in the comoving frame of the plasma fluid, 

€ = p-pB: pB=nBmc^ , (19) 

and using the law of baryon-number conservation, from Eq. (|18fl we have 

{eU''),,+p{U''),,=Q . (20) 

Recalhng that ^ = y([/'").^, where V is the comoving volume and t is the proper 
time for the plasma fluid, we have along each flow line 

dE 
d7 



d{Ve) ^ dy_ 
dr ^ dr 



n 



(21) 



where E = Ve is the total proper internal energy of the plasma fluid. We express 
the equation of state by introducing a thermal index T{p, T) 

r = 1 + - . (22) 

e 

We now turn to the second set of governing equations describing the evolution 
of the e^e^ pairs. Letting n,.- and ne+ be the proper number densities of electrons 
and positrons associated with pairs and the proper number densities of ionized 
electrons, we clearly have 



''pair 7 



= Zn 



B, 



(23) 
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where ripair is the number of e+e^ pairs and Z the average atomic number i < 
Z < 1 (Z = 1 for hydrogen atom and Z — ^ for general baryonic matter). The rate 
equation for electrons and positrons gives, 



1 



= aU[K-(r)+n^_(T))n,+ (T) 
-K-+n^_K+], (24) 

K-C/^);,, - K-C/*),t + ^{r^n,-Un,r 

^av[n^-{T)n^+{T)-n^-n^+], (25) 

^mj[nl^{T)n,+ {T)-nl,n,+ ], (26) 

where av is the mean of the product of the annihilation cross-section and the 
thermal velocity of the electrons and positrons, ng± (T) are the proper number 
densities of electrons and positrons associated with the pairs, given by appropriate 
Fermi integrals with zero chemical potential, and ri^_ (T) is the proper number 
density of ionized electrons, given by appropriate Fermi integrals with non-zero 
chemical potential fie at an appropriate equilibrium temperature T. These rate 
equations can be reduced to 

(n,±[/'^).^ = (ne±[/*),t + ^(r^n.i [/'■),, 

= CTu[ne-(T)ne+(T) - rie-ne+j, (27) 
{nlU^;^ = {nl-U'),^ + l{r^nlun,r = 0, (28) 

Frac^^ = ^^. (29) 

Equation H28|) is just the baryon-number conservation law and (|29|l is a rela- 
tionship satisfied by ne± , ng± (T) and n^. , n^. (T). 

The equilibrium temperature T is determined by the thermalization processes 
occurring in the expanding plasma fluid with a total proper energy density p gov- 
erned by the hydrodynamical equations pt)ll7ll8|l . We have 

P = Pf + Pe+ + Pe- + pI- + PB, (30) 

where is the photon energy density, pB — raBC^riB is the baryonic mass den- 
sity which is considered to be nonrelativistic in the range of temperature T under 
consideration, and pg± is the proper energy density of electrons and positrons pairs 
given by 

Pe± = "'''Lx Pe± {T), (31) 
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where n^± is obtained by integration of Eq. I|27|l and (T) is the proper energy 
density of electrons(positrons) obtained from zero chemical potential Fermi integrals 
at the equilibrium temperature T. On the other hand Pg_ is the energy density of 
the ionized electrons coming from the ionization of baryonic matter 

Pe-=;^Pe-(T^)' (32) 

where is obtained by integration of Ea. (|28|l and p^-{T) is the proper energy 
density of ionized electrons obtained from an appropriate Fermi integral of non-zero 
chemical potential /ig at the equilibrium temperature T . 

Having intrinsically defined the equilibrium temperature T in Eq.JSOI), we can 
also analogously evaluate the total pressure 

P^P-l +Pe+ +Pe- +P\- +PB, (33) 

where p^ is the photon pressure, Pe± ^-nd p\- are given by 

Pe± = ^^^^^Pe±(T), (34) 
ne± (T) 

the pressures Pe± {T) are determined by zero chemical potential Fermi integrals, 
and pI. (T) is the pressure of the ionized electrons, evaluated by an appropriate 
Fermi integral of non-zero chemical potential /Xe at the equilibrium temperature 
T. In Ea. (|33|) . the ion pressure pb is negligible by comparison with the pressures 
g± g- (T), since baryons and ions are expected to be nonrelativistic in the range 
of temperature T under consideration. Finally using Eas. H3UI33|l we compute the 
thermal factor F of the equation of state (|22|l . 

It is clear that the entire set of equations considered above, namely Eas. H16ll7ll8|) 
with equation of state given by Ea. (|^ and the rate equation (|?7|l . have to be in- 
tegrated satisfying the total energy conservation for the system. The boundary 
conditions adopted here are simply purely ingoing conditions at the horizon and 
purely outgoing conditions at radial infinity. The calculation is initiated by deposit- 
ing a proper energy density between the Reissner-Nordstrom horizon radius 
r-|_ and the dyadosphere radius r^s, following the approximation presented in Fig El 
The total energy deposited is given by Ea. fT^ . 



4. The equations leading to the relative space-time 
transformations 

In order to relate the above hydrodynamic and pair equations with the observations 
we need the governing equations relating the comoving time to the laboratory time 
corresponding to an inertial reference frame in which the EMBH is at rest and 
finally to the time measured at the detector, which must also include the effect 
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of the cosmological expansion. These transformations have been the object of the 
relative space-time transformations (RSTT) ParadigmP 

For signals emitted by a pulse moving with velocity v in the laboratory frameP 
we have the following relation between the interval of arrival time Ai^ and the 
corresponding interval of laboratory time At (see Fig. EJ: 

Ma={to + M+^^^-^^-{to + ^^^At~'-. (36) 

For simplicity in what follows we indicate by the interval of arrival time mea- 
sured from the reception of a light signal emitted at the onset of the gravitational 
collapse. Analogously, t indicates the laboratory time interval measured from the 
time of the gravitational collapse. In this case, Ea. H36|l can be written simply as: 

t.^t-^-^t-^>^nj^^^^ (37) 

c c 

where the dyadosphere radius rds is the value of r at i = 0. We consider here only 
the photons emitted along the line of sight from the external surface of the pulse. 
The arrival time spreading due to the angular de penden ce and that due to the 
thickness of the pulse will be considered elsewhere! I I The solution of Eq. (|37|l 
has the expansion: 

ta=t-—t-\—e~..., (38) 

so the relation between ta and t is in general highly nonlinear. 

If and only if the expansion of the pulse is such that r [t) ~ vt with w ~ c, Eq. (|37|l 
can be written, neglecting rds, in the following simplified form (see Fig. 110(1 : 

V c) + ^ ' 

This formula has been uncritically and widely applied in all articles dealing with 
GRBs. It is clear, however, that the knowledge of ta, which is indeed essential for 
any physical interpretation of GRB data, depends on the definite integral given 
in Ea. (|37|l whose integration limits in the laboratory time extend from the onset 
of the gravitational collapse to the time t relevant for the observations. Such an 
integral is not generally expressible as a simple linear relation or even by any explicit 
analytic relation since we are dealing with processes with variable gamma factor 
unprecedented in the entire realm of physics (see Figs. Inland Fig. I1U|) . Any lin ear 
approximation of the kind given in Eq. 139() with 7 constant or changing with timi^ 
misses a crucial feature of the GRB process and is therefore erroneous in this 
context. 

To relate the time in the laboratory frame to the time in the detector frame 
we have to do one additional step: the two frames are related by a transformation 
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which is a function of the cosmological expansion. We recall that the geometry of 
the space-time of the universe is described by the Robertson- Walker metric: 

ds^ = - n^(t) — -t- r^d^^ + r^sin-d^dip^ , (40) 

\ 1 — fcr^ J 

where TZ (t) is the cosmic scale factor and fc is a constant related to the curvature 
of the three-dimensional space {k = 0, -fl,— 1 corresponds to flat, close and open 
space respectively). The wavelength of an electromagnetic wave travelling from the 
point Pi(ti, ri, (/9i) to the point Po{to,ro,^9o, ifio) where the observer is located 
is related to the red-shift parameter z by 

z = , (41) 

M 

where Ao is the wavelength of the radiation for the observer and Ai for the emitter. 
We have the following general relation: 

l + z = (l + z0(l + zo)(l + z,), (42) 

where z is the total redshift due to the motion of the source Zg, the motion of 
the observer Zg and the cosmological redshift z^. In the following we will assume 
Zo << 1 and << 1 so z = z„. In terms of the scale factor TZ (t) the relation (|41|l 
gives 

^ = lM = l + - = - (43) 

where uji and wq £^re the frequencies associated to Ai and Aq respectively. This 
frequency ratio then relates the time elapsing at the source with the time elapsing 
at the detector due to the cosmological expansion. 

We can now define the corrected arrival time fj^ measured at the detector, which 
is related to ta by 

tt^t,{l + z), (44) 

where z is the cosmological redshift of the GRB source. In the case of GRB 991216 
we have z ~ 1.00. 

The observed flux is the flux which crosses the surface 47r(7?. {to) r)^ but this 
flux is lower by a factor 1 -\- z due to the redshift energy of the photons and by 
another factor 1 + z due to the fact that the number of photons at reception is less 
than the number at emission. Thus we can define a luminosity distance by: 

dl^Tiy{l + zf. (45) 



Then the observed flux is related to the absolute luminosity of the GRB by the 
following relation: 
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where the luminosity distance cIl is simply related to the proper distance dp — TZqI' 
hy (11 — dp{l + z) . The observed total fluence / is related to the total energy E of 
the GRB by the following relation: 

/ = + (47) 

Then the cosmological effect is taken into account by the definition of the proper 
distance TZof which depends on the cosmological parameters: the Hubble constant 
Ho = TZ (to) /"R- (to) at time to and the matter density po or Qm ~ Po/ Pcrit, where 

Pcrit — g^Q- 

The computation of the proper distance is then simply given by the relation : 



c 



where F{z) = y/flMil + z)^. 

In the case of the Friedman flat universe, Qm ~ 1 and we have: 

1 



tin 



1 



(49) 



So the measurement of the redshift gives us the luminosity distance via a cos- 
mological scenario. With the measurement of the flux we can deduce the proper 
luminosity of the burst and from the measurement of the total fluence the total 
energy so we are then able to find the Edya- 

5. The numerical integration of the hydrodynamics and the rate 
equations 

5.1. The Livermore code 

A computer codJ^^EHI has been used to evolve the spherically symmetric general 
relativistic hydrodynamic equations starting from the dyadosphere!^ 

We define the generalized gamma factor 7 and the radial 3-velocity in the lab- 
oratory frame 

jjr 



From Eas.l(5l we then have 

1 1 

gtt 

Following Ea. H19|l . we also define 

E= €^j, D = pBl, and p= pj (52) 
so that the conservation law of baryon number 1)16(1 can then be written as 

dD a d 

^ = --^(-W). (53) 
at r-^ or a 
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Ea. H18|) then takes the form, 



dE_ 
'dt 



a d ,r 



dr 



d-f 
'dt 



a d 

J.2 Qj, 



(54) 



Defining the radial momentum density in the laboratory frame 

Sr = a{p + p)U^Ur^{D + TE)Ur, (55) 

we can express the radial component of the energy-momentum conservation law 
given in Eq.lQTJ by 



dSr 



or a 



a- 



dp 



a . , 



d9u . 
or 



dr 

dr 



or a Or 



M 



D + TE 

7 



r\2-\ 



(56) 



In order to determine the number-density of e^e^ pairs, we turn to Ea. (|27|l . 
Defining the e"'"e~-pair density in the laboratory frame N^± = ^n^± and N^± (T) = 
7ne±(T), where the equilibrium temperature T has been obtained from Eas. H30|l 
and H^lfl . and using Ea. (|51|l . we rewrite the rate equation given by Ea. pTfl in the 
form 

dN^± a d ,r'^ _ ^ , ■? 

(57) 



dt 



These equations are integrated starting from the dyadosphere distributions given 
in Fig. 1171 and assuming as usual ingoing boundary conditions on the horizon of the 
EMBH. 



5.2. The Rome code 

In the following we recall a z erot h order approximation of the fully relativistic equa- 
tions of the previous sectionS^ 

(i) Since we are mainly interested in the expansion of the e"*"e~ plasma away from 
the EMBH, we neglect the gravitational interaction. 

(ii) We describe the expanding plasma by a special relativistic set of equations. 

(iii) In contrast with the previous treatment where the evolution of the density pro- 
files given in Fig. El are followed in their temporal evolution leading to a pulse-like 
structure, selected geometries of the pulse are a priori adopted and the correct one 
validated by the complete integration of the equations given by the Livermore codes. 

In analogy to Ea. H21|l . from Eg. ljlGII we have along each flow line in the general 
case in which baryonic matter is present 

%^=0. (58) 
dr 



On the structure of the burst and afterglow of Gamma-Ray Bursts I: the radial approximation 37 



For the expansion of a shell from its initial volume AVo to the volume Al/^, we 
obtain 

n| ^ AT/ ^ AV7(r) 

UB AK AVo7o(r)' 
where AV is the volume of the shell in the laboratory frame, related to the proper 
volume AV^ in the comoving frame by AV^ = 7(r)AV, where 7(r) defined in Ea. (|50|l 
is the gamma factor of the shell at the radius r. 

Similarly from Ea. (|21|l . using the equation of state (|22ll . along the flow lines we 
obtain 

dlne + rdlnV" = 0. (60) 
Correspondingly we obtain for the internal energy density e along the flow lines 



r 

Co 



(61) 



Ay y _ ( ^\ ( li-r) 

where the thermal index T given by H22|l is a slowly-varying function with values 
around 4/3. It can be computed for each value of e,p as a function of AF. 

The overall energy conservation requires that the change of the internal proper 
energy of a shell is compensated by a change in its bulk kinetic energy. We then 
havJSl 

dK =[-i{r)~l]{dE + pBdV). (62) 

In order to model the relativistic expansion of the plasma fluid, we assume that 
E and D as defined by Ea. (|52|l are constant in space over the volume AF. As a 
consequence the total energy conservation for the shell implies— 

(eo + p°b)iI{t)AVo = {e + pb)i^{t)/\V, (63) 
which leads the solution 



V (e + ps)AV 

Corresponding to Ea. H57|l we obtain the equation for the evolution of the 
number-density as seen by an observer in the laboratory frame 

|(A^..) = -^e.^^ + -7-i-^(iV|4T) ~ Nl,) . (65) 

Eas. H59|l . H61|l . (|64l) and l|65() are a complete set of equations describing the rel- 
ativistic expansion of the shell. If we now turn from a single shell to a finite 
distribution of shells, we can introduce the average values of the proper internal- 
energy, baryon-mass, baryon- number and pair- number densities (e, pB, nB, n^±) and 
E = ^e, D = jpb, = 7(r)ne± for the PEM-pulse, where the average 7-factor 
is defined by 

l=\j^l(r)dV, (66) 
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and V is the total volume of the shell in the laboratory frame. The corresponding 
equations are given in Ruffini, Salmonson, Wilson & Xue (f 999)P^Having defined 
all its governing equations we can now return to the description of the different eras 
of the GRB phenomena. 



6. The era I: the PEM pulse 

We have assumed that, following the gravitational collapse process, a region of very 
low baryonic contamination exists in the dyadosphere all the way to the remnant 
of the progenitor star. 

Recalling Eq.ljS)) the limit on such baryonic contamination, where pB^ is the 
mass-energy density of baryons, is given by 



< mprie+e- W = 3.2 • 10^ 



{g I cm? 



Near the horizon 



this gives 



PB, < TOpne+e- (r) = 1.86 • lO^"* ( - ) (.g/cm^) . 



and near the radius of the dyadosphere Vds 



PBc mpTie+e- {r) — 3.2 • 10^ 



1-1 — 

Tds 



{gjcw? 



(67) 



(68) 



(69) 



Such conditions can be easily satisfied in the collapse to an EMBH, but not neces- 
sarily in a collapse to a neutron star. 

Consequently we have solved the equations governing a plasma composed solely 
of e^e^-pairs and electromagnetic radiation, starting at time zero from the dyado- 
sphere configurations corresponding to constant density in Fig. 1171 The Livermore 
codeP 

has shown very clearly the self organization of the expanding plasma in a 
very sharp pulse which we have defined as the pair-electromagnetic pulse (PEM 
pulse), in analogy with the EM pulse observed in nuclear explosions. In order to 
further examine the structure of the PEM pulse with the simpler procedures of 
the Rome codes we have assume three alternative patterns of expansion of the 
PEM pulse on which to try the simplified special relativistic treatment and then 
compared the results with the fully general relativistic hydrodynamical results: 

• Spherical model: we assume the radial component of the four- velocity Ur{r) = 
U^, where U is the radial component of the four- velocity at the moving outer 
surface r = Ti{t) of the PEM pulse and the 7-factor and the velocity Vr are 

3 



7 



2[/(n-c/^)2 -uii + u^y- 



ln(f/ -t- ^/T+TP) 



7 



(70) 
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this distribution expands keeping an uniform density profile wliicli decreases 
with time similar to a portion of a Friedmann Universe. 

• Slab 1: we assume U{r) — Ur = const., the constant width of the expanding 
slab V = Ro in the laboratory frame of the PEM pulse, while 7 and Vr are 

7 = v/TTW, K^^; (71) 
7 

this distribution does not need any averaging process. 

• Slab 2: we assume a constant width R2 — Ri = Ro of the expanding slab in the 
comoving frame of the PEM pulse, while 7 and Vr are 

7= Vl + C/2(f), K = (72) 

7 

This distribution needs an averaging procedure and Ri < r < R2, i.e. f is an 
intermediate radius in the slab. 

These different assumptions lead to three different distinct slopes for the mono- 
tonically increasing 7-factor as a function of the radius (or time) in the laboratory 
frame, having assumed for the energy of dyadosphere E'dya — 3.1 x 10^^ erg (see 
Fig. I18|) . In principle, we could have an infinite number of models by defining ar- 
bitrarily the geometry of the expanding fluid in the special relativistic treatment 
given above. To find out which expanding pattern of PEM pulses is the physically 
realistic one, we need to compare and contrast the results of our simplified mod- 
els (performed in Rome) with the numerical results based on the hydrodynamic 
Eas. l|53l54l56|) (obtained at Livermorc).^^ Details of the iterative method used to 
solve the special relativistic equation can be found in Ruffini, Salmonson, Wilson 
& Xue (1999).42! 

It is manifest from the results (see Fig. II8II that the slab 1 approximation (con- 
stant thickness in the laboratory frame) is in excellent agreement with the Livermore 
results (open squares). 

The remarkable validation of the special relativistic treatment of the PEM 
pulse,t-^ allows us to easily estimate the related quantities of physical and astro- 
physical interest in the model, like the e+e~-pair densities as a function of the 
laboratory time, the temperature of the plasma in the comoving and laboratory 
frames, the reheating ratio as a function of the e+e~-pair annihilation for a variety 
of initial conditions 

7. The era II: the interaction of the PEM pulse with the remnant 
of the progenitor star 

The PEM pulse expands initially in a region of very low baryonic contamination 
created by the process of gravitational collapse. As it moves further out the baryonic 
remnant (see Fig. ^ of the progenitor star is encountered. As discussed in section 
1^ below, the existence of such a remnant is necessary in order to guarantee the 
overall charge neutrality of the system: the collapsing core has the opposite charge 
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Fig. 18. Gamma factor as a function of radius. Three models for the expansion pattern of the 
PEM-pulse are compared with the results of the one dimensional hydrodynamic code for an energy 
of dyadosphere E^y^ = 3.1 X 10^^ erg. The 1-D code has an expansion pattern that strongly 
resembles that of a shell with constant thickness in the laboratory frame. 

of the remnant and the system as a whole is clearly neutral. The number of extra 
charges in the baryon ic remnant negligibly affects the overall charge neutrality of 
the PEM pulse.CSESI 

The baryonic matter remnant is assumed to be distributed well outside the 
dyadosphere in a shell of thickness A between an inner radius rjn and an outer 
radius rout = Tin + A at a distance from the EMBH at which the original PEM 
pulse expanding in vacuum has not yet reached transparency. For the sake of an 
example we choose 

Tin = lOOrds, A = lOrd.. (73) 
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The total baryonic mass Mb = Nsmp is assumed to be a fraction of the dyado- 
sphere initial total energy (-Edya)- The total baryon- number Nb is then expressed 
as a function of the dimensionless parameter B given by 

S = ^^, (74) 

where i? is a parameter in the range 10~^ — 10~^ and nip is the proton mass. 
We shall see below the paramount importance of B in the determination of the 
features of the GRBs. We will see in section |51 the sense in which B and E^ya can 
be considered to be the only two free parameters of the EMBH theory for the entire 
GRB family, the so called "long bursts" . We shall see in section ^2 that for the so 
called "short bursts" the EMBH theory depends on the two other parameters ft, 
^, since in that case B = 0. The baryon number density n°g is assumed to be a 
constant 

- o -o -Q 2 

ng^—, Pb = rnpTiBC . (75) 
Vb 

As the PEM pulse reaches the region ri„ < r < rout, it interacts with the 
baryonic matter which is assumed to be at rest. In our simplified quasi-analytic 
model we make the following assumptions to describe this interaction: 

• the PEM pulse does not change its geometry during the interaction; 

• the collision between the PEM pulse and the baryonic matter is assumed to be 
inelastic, 

• the baryonic matter reaches thermal equilibrium with the photons and pairs of 
the PEM pulse. 

These assumptions are valid if: (i) the total energy of the PEM pulse is much 
larger than the total mass-energy of baryonic matter Mb, 10^^ < B < 10^^, (ii) 
the ratio of the comoving number density of pairs and baryons at the moment of 
collision Ue+e-Ms is very high (e.g., 10^ < n^+e- /n°g < 10^^) and (iii) the PEM 
pulse has a large value of the gamma factor (100 < 7). 

In the collision between the PEM pulse and the baryonic matter at rout > r > r^n 
, we impose total conservation of energy and momentum. We consider the collision 
process between two radii r2,ri satisfying Tout > r2 > ri > r-^ and r2 — ri ^ A. 
The amount of baryonic mass acquired by the PEM pulse is 

AM = ^f(r3-r?), (76) 

where Mb/Vb is the mean-density of baryonic matter at rest. The conservation of 
total energy leads to the estimate of the corresponding quantities before (with "o" ) 
and after such a collision 

(Fe-o + p^)7o Vo + AM = (Fe + + ^ + TAl)fV, (77) 
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where Ae is the corresponding increase of internal energy due to the coUision. 
Similarly the momentum-conservation gives 

(reo + P°b)1oU:V, = (re + PS + ^ + rAe)7C/,V, (78) 



where the radial component of the four-velocity of the PEM pulse is U° — 
and r is the thermal index. We then find 



r 



, -o, 7ot/,°Vo AAf , 



(79) 



a _ 7o AA/ 

^"7^^' ^ (re-o + p^)7oC/,°Vo' ^^^^ 

These equations determine the gamma factor 7 and the internal energy density 
Z—Zo + Ae in the capture process of baryonic matter by the PEM pulse. 

The effect of the collision of the PEM pulse with the remnant leads to the 
following r esultPl as a function of the B parameter defined in Ea. H74|l : 

1) an abrupt decrease of the gamma factor given by 

where 70 is the gamma factor of the PEM pulse prior to the collision and B is given 
by Eq.lf7l|. 

2) an increase of the internal energy in the comoving frame EcoU developed in the 
collision given by 



Edya 7o \T 

3) a corresponding reheating of the plasma in the comoving frame but not in the 
laboratory frame, an increase of the number of e+e^ pairs and correspondingly an 
overall increase of the opacity of the pulse. See details in section [TUl 



8. The era III: the PEMB pulse 

After the engulfment of the baryonic matter of the remnant the plasma formed 
of e+e~-pairs, electromagnetic radiation and baryonic matter expands again as a 
sharp pulse, namely the PEMB pulse. The calculation is continued as the plasma 
fluid expands, cools and the e"'"e~ pairs recombine until it becomes optically thin: 

/ dr(7ie± +ZnB)o-T (83) 

JR 

where ctt — 0.665 • 10^^'^cm^ is the Thomson cross-section and the integration 
is over the radial interval of the PEMB pulse in the comoving frame. We have 
first explored the general problem of the PEMB pulse evolution by integrating 
the general relativistic hydrodynamical equations with the Livermore codes, for a 



On the structure of the burst and afterglow of Gamma-Ray Bursts I: the radial approximation 43 



total energy in the dyadosphere of 3.1 x 10^'* erg and a baryonic shell of thickness 
A = lOrds at rest at a radius of lOOrds and S ~ 1.3 • 10"^. 

In total analogy with the special relativistic treatment for the PEM pulse, pre- 
sented in section (see also Ruffini, Salmonson, Wilson & Xue, 19954%, we obtain 
for the adiabatic expansion of the PEMB pulse in the constant-slab approximation 
described by the Rome codes the following hydrodynamical equations with ps ^ 





V 


TLB 


Vo ~ 


eo _ 




e 




7 = 





V7 

Vo7o ' 
r 



{TI+pb)V ' 



(84) 
(85) 

(86) 
(87) 



In these equations (r > rout) the comoving baryonic mass- and number densities 
are pB — Mb/V and tIb — Nb/V, where V is the comoving volume of the PEMB 
pulse. 

We compare and contrast (see Fig. (S)) the bulk gamma factor as computed from 
the Rome and Livermore codes, where excellent agreement has been found. This 
validates the constant-thickness approximation in the case of the PEMB pulse as 
well. On this basis we easily estimate a variety of physical quantities for an entire 
range of values of B. 

For the same EMBH we have considered five different cases: a shell of baryonic 
mass with (1) B ~ 1.3 • lO""; (2) B ~ 3.8 • IQ-*; (3) B ~ 1.3 • IQ-^; (4) B ~ 
3.8 • 10^"^; (5) B ~ 1.3 • 10"'^). The results of the integration given in detail in 
Ruffini, Salmonson, Wilson & Xue (2000)^'^ show that for the first parameter range 
the PEMB pulse propagates as a sharp pulse of constant thickness in the laboratory 
frame, but already for i? ~ 1.3 • 10~^ the expansion of the PEMB pulse becomes 
much more complex and the constant-thick ness approximation ceases to be valid; 
see Ruffini, Salmonson, Wilson & Xue (2000 for details. 

It is particularly interesting to evaluate the final value of the gamma factor of 
the PEMB pulse when the transparency condition given by Ea. l83|l is reached as 
a function of B, see Fig. 1191 For a given EMBH, there is a maximum value of the 
gamma factor at transparency. By further increasing the value of B the entire E^ya 
is transferred into the kinetic energy of the baryons; see also section [TTl Details are 
given in Ruffini, Salmonson, Wilson & Xue (2000). 

In Fig. |2D1 we plot the gamma factor of the PEMB pulse versus the radius for 
different amounts of baryonic matter. The diagram extends to values of the radial 
coordinate at which the transparency condition given by Eq.lIH3 is reached. The 
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Fig. 19. The gamma factor (the sohd line) at the transparent point is plotted as a function of 
the B parameter. The asymptotic value (the dashed line) E^y^/{AlBC^) is also plotted. 



"asymptotic" gamma factor 



Tasym 



E, 



dya 



(88) 



is also shown for each curve. The closer the gamma value approaches the "asymp- 
totic" value H88|l at transparency, the smaller the intensity of the radiation emitted 
in the burst and the larger the amount of kinetic energy left in the baryonic matter. 



9. The identification of the free parameters of the EMBH theory 

Within the approximation presented in section|21the EMBH is characterized by two 
parameters: /i and ^. The energy of the dyadosphere is expressed in terms of these 
two parameters by Ea. l(T^ . 

There is an entire family of EMBH solutions with different values of /i and ^ 
corresponding to the same value of Edya (see Fig. ll6|l . These solutions are physically 
different with respect to the density of electron-positron pair distributions given by 
Eq.Q, as well as to their energy density given by Eq.lQTJ. A clear example of 
such a degeneracy is given in Fig. 1211 where the two limiting energy density profiles 
approximating the dyadosphere as introduced in Fig. 1 171 are given for three different 
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Fig. 20. The gamma factors are given as functions of the radius in units of the dyadosphere 
radius for selected values of B for the typical case Edyo. = 3.1 X 10^* erg. The asymptotic values 
7agy^ = E^y^l{MBC^) = 10*, 10^, 10^ are also plotted. The collision of the PEM pulse with the 
baryonic remnant occurs at rfr^^ = 100 where the jump occurs and the PEMB pulse starts. 



EMBH configurations corresponding to the same value of E'dj/a = 3.1 x 10 erg. 
The three configurations correspond respectively to the three different pairs (/i,^): 
(10,0.76), (102,0.27), (10^,0.10). 

The corresponding dynamical evolution of the P EM pulse introduced in sec- 
tion Eland RufRni, Salmonson, Wilson & Xue (1999 is clearly different in the 
three cases. It is remarkable that when the collision with the remnant of the pro- 
genitor star is considered all these differences disappear. As usual (see section^ we 
describe the baryonic content of the remnant by the parameter B. The PEMB pulse 
generated after the collision with the baryonic matter depends uniquely on the two 
parameters E^ya and B. In Fig. l22l the temperature in the laboratory frame is given 
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Fig. 21. Three different dyadosplieres corresponding to the same value of E^iya = 3.1 X 10 erg 
and with different values of the two parameters and ^ are given. The three different configurations 
are markedly different in their spatial extent as well as in their energy-density distribution. 



for the PEM pulse and the PEMB pulse corresponding to the three configurations 
of Figl^and B = Ax 10~^. It is clear that while for the PEM pulse era the three 
configurations are markedly different, they do converge to a common behaviour in 
the PEMB pulse era. 

If we turn now to the effect of the distance between the EMBH and the baryonic 
remnant, we see that this degeneracy is further extended: while the three PEM 
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Fig. 22. The temperature of the plasma during the PEM pulse and PEMB pulse eras, measured 
in the laboratory frame, corresponding to the three configurations presented in Fig. 1211 is given 
as a function of the laboratory time. The three different curves converge to a common one in the 
PEMB pulse era, which is therefore only a function of the Edya and B. The difference among 
the three curves in the early part of the PEMB pulse follows from having located the baryonic 
matter at a distance of 50(rds — which is different in the three cases. Such difference become 
negligible at large distances in the later phases of the evolution. 



pulse eras are quite different, the common PEMB pulse era is largely insensitive 
to the location of the baryonic remnant, see Fig. 123 We have plotted the three 
gamma factors in the PEM pulse era corresponding to the different configurations 
of Fig. 1211 and B = 10^^, in the two cases the baryonic remnant is positioned at 
different distances from the EMBH. 

If the PEM pulse has reached extreme relativistic regimes, the common value 
7co;/ to which the three gamma factors drop in the collision with the baryonic matter 
of the remnant can be simply expressed by the large gamma limit of Eq. (|81ll 

while the internal energy Ecoii developed in that collision is simply given by the 
corresponding limit of Eq. H82(l 



coll 



^dya 



-B + \/B^ + 2B . 



(90) 



This approximation applies when the final gamma factor at the end of the PEM 
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pulse era is larger than ^coU , upper panel in Fig. |2S1 
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Fig. 23. The gamma factors for the three configurations considered in Fig. 1211 are given as a 
function of the radial coordinate in the laboratory frame. The two figures correspond to a baryonic 
remnant positioned respectively at ri„ = 50(r;jg — r+) (above) and at ri„ = 5{r^g — *"+)• Again 
the convergence to a common behaviour, uniquely a function of E^iya and B for the late stages of 
the PEMB pulse, is manifest. 



Turning from these general considerations to the GRB data, this degeneracy in 
the PEMB pulse eras and their dependence on only two parameters Edya and B has 
far reaching astrophysical implications for the identification of the source of GRBs. 
As we will see in the conclusions all the information obtainable from GRBs with 
a large value of the parameter B will lead to the determination of the above two 
parameters. An entire family of degenerate astrophysical solutions in the range of 
charges and masses given in Fig. E| are possible. The direct knowledge of the mass 
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and charge of the EMBH can only be gained from the PEM pulse or from GRBs 
with very small values oi B — the so called "short bursts", see section ITTl and the 
conclusions. 



10. The approach to transparency: the thermodynamical 
quantities 

As the condition of transparency expressed by Eq. 1)83(1 is reached the injector phase 
terminates. The electromagnetic energy of the PEMB pulse is released in the form 
of free-streaming photons — the proper GRB. The remaining energy of the PEMB 
pulse is released as an accelerated-baryonic-matter (ABM) pulse. 

We now proceed to the analysis of the approach to the transparency condition. 
It is then necessary to turn from the pure dynamical description of the PEMB 
pulse described in the previous sections to the relevant thermodynamic parameters. 
Also such a description at the time of transparency needs the knowledge of the 
thermodynamical parameters in all previous eras of the GRB. 

As above we shall consider as a typical case an EMBH of Edya = 3.1 x 10^'' 
erg and B = 10~^. The considerations will refer to a dyadosphere configuration 
described by the two limiting approximations shown in Fig. 1171 

One of the key thermodynamical parameters is represented by the temperature 
of the PEM and PEMB pulses. It is given as a function of the radius both in the 
comoving and in the laboratory frames in Fig. |^ Before the collision the PEM 
pulse expands keeping its temperature in the laboratory frame constant while its 
temperature in the comoving frame falls .1^. In fact Eas. H63l64|l are equivalent to 

where the baryon mass-density is pB =0 and the thermal energy-density of photons 
and e~''e~-pairs is e = (T_bT^(1 + fe+e-), <^b is the Boltzmann constant and fe+e- 
is the Fermi-integral for e"*" and e~. This leads to 

e-f^V = Edy^, T^"/^V = const. (92) 

Since e"*" and in the PEM pulse are extremely relativistic, we have the equation 
of state p ~ e/3 and the thermal index H22|l F ~ 4/3 in the evolution of PEM pulse. 
Ea. ((92|l is thus equivalent to 

T^jV ~ const. (93) 

These two equations ()91(l and ((93|l result in the constancy of the laboratory tem- 
perature in the evolution of the PEM pulse. 

It is interesting to note that Eas. ((92|l and ((93(1 hold as well in the cross-over 
region where T ~ mgC^ and e^e^ annihilation takes place. In fact from the conser- 
vation of entropy it follows that asymptotically we have 
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exactly for the same reasons and physics scenario discussed in the cosmological 
framework by Weinberg, see e.g. Eq. (15.6.37) of Weinberg (1972). The same con- 
siderations when repeated for the conservation of the total energy ejV = e7^V 
following from Eq. I|9HI then lead to 

{VTS)T<m.c^ _ 11 

The ratio of these last two quantities gives asymptotically 

where To is the initial average temperature of the dyadosphere at rest. 

During the collision of the PEM pulse with the remnant we have an increase in 
the number density of e~^e~ pairs (see Fig. I25II . This transition corresponds to an 
increase of the temperature in the comoving frame and a decrease of the temperature 
in the laboratory frame as a direct effect of the dropping of the gamma factor (see 
Fig.EOI). 

After the collision we have the further acceleration of the PEMB pulse (see 
Fig. I20|l . The temperature now decreases both in the laboratory and the comoving 
frame (see Fig. I24() . Before the collision the total energy of the e+e~ pairs and the 
photons is constant and equal to i?dya- After the collision 

+ -Bg+e- + -Ephotons, (97) 

which includes both the total energy -Eg+e- + -Ephotons of the nonbaryonic compo- 
nents and the kinetic energy i^Baryons of the baryonic matter 

EBaryon. =PSV^(7-1)- (98) 

In Fig. 1201 we plot both the total energy Ee+e- + Ephotons of the nonbaryonic 
components and the kinetic energy i?Baryons of the baryonic matter as functions of 
the radius for the typical case Edya = 3.1 X lO^'' erg and B = 10"^ Further details 
are given in Ruffini, Salmonson, Wilson & Xue (2000"^^ 



11. The P-GRBs and the "short bursts". The end of the injector 
phase. 

We now analyze the approach to the transparency condition given by Eq. H83|) . For 
selected values of B we give the energy Ep-qrb of the P-GRB, and i?Baryons of the 
ABM pulse. We clearly have 

Edya — Ep-GRB + EBaryons • (99) 

Taking into account the results shown in Figs. I24II26I we can repeat all the 
considerations for selected values of B. We shall examine values of B ranging from 
B = 10^^ only up to _B = 10^^: for larger values of B our constant slab approxi- 
mation breaks down. We will see in the following that this range does indeed cover 
the most relevant observational features of the GRBs. 
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Fig. 24. The temperature of the plasma in the comoving frame T'(MeV) (the solid line) and in 
the laboratory frame 7T' (the dashed line) are plotted as functions of the radius in the unit of 
the dyadosphere radius r^g. 
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Fig. 25. The number densities n^+g- (T) (the solid line) computed by the Fermi integral and 
rig+g- (the dashed line) computed by the rate equation (see section|2l are plotted as functions of 
the radius. T' <C mf,c^ , two curves strongly divergent due to e+e~-pairs frozen out of the thermal 
equilibrium. The peak at r ~ lOOrds is due to the internal energy developed in the collision. 



As clearly shown in Fig.l^Ulbotli the final value of the gamma factor and the ra- 
dial coordinate at which the transparency condition is reached depend very strongly 
on B. Therefore a strong dependence on B is also found in the relative values of 

Ep-GRB and i^Baryons- 

We are now finally ready to give in Fig. |^ the crucial diagram representing the 
values of Ep-qi^b and -EBaryons in units of the E^ya as functions of B. This diagram, 
a universal one, is very important and is essential for the understanding of the GRB 
structure. 




Fig. 26. The energy of the non bciryonic components of the PEMB pulse (the solid line) and the 
kinetic energy of the baryonic matter (the dashed line) in unit of the total energy are plotted as 
functions of the radius in the unit of the dyadosphere radius r^s ■ 
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We find that for small values of B (around 10 ^) almost all the Edya is emitted 
in the P-GRB (see also our previous paper Ruffini, Salmonson, Wilson & Xue, 
and very little energy is left in the baryons. While for i? ~ 10 roughly 
only 10~^ of the total initial energy of the dyadosphere is radiated away in the 
P-GRB and almost all energy is transferred to the baryons. 

This behaviour is at the heart of the fundamental difference between the so 
called short hursts and long bursts. We have proposec0 that the short bursts must 
be identified with the P-GRBs in the case of very small B. There are a variety of 
reasons supporting this identification: 



(1) For small values of B, i?Baryons is negligible, see Fig. and consequently the 
intensity of the afterglow is also negligible and the entire energy E^ya is re- 
leased into the P-GRB. This is clearly consistent with the absence of observed 
afterglows in the short bursts. 

(2) The temperature of the P-GRB in the laboratory frame jT at the transparency 
point is a strongly decreasing function of B, see Fig. [T] 7T is related to the 
energy corresponding to the peak of the ph oton -number spectrum, as described 
in Ruffini, Salmonson, Wilson & Xue This is also in very good agree- 
ment with the observed decrease of the hardness ratio between the short bursts 
and the long bursts^^ 

(3) The time Tqq, the duration of 90% of the energy emission as used in the cur- 
rent literature and discussed in Ruffini, Salmonson, Wilson & Xue (2000^P^ is 
plotted in Fig. [2Z|for selected values of Edya and for different values of B. 



Before concluding a word of caution is needed about how to use the above results: 
all these considerations are based on the drastic approximations in the description 
of the dyadosphere presented in section [51 see also Fig. [21 This treatment is very 
appropriate in estimating the general dependence of the energy of the P-GRB, the 
kinetic energy of the ABM pulse and consequently the intensity of the afterglow. 
Especially powerful is the establishment of the dependence of Ep-grb and i^Baryons 
on B (see Fig.[^. As we will see in the next sections, this approximation is similarly 
powerful in determining the overall time structure of the GRB and especially the 
time of the release of the P-GRB with respect to the moment of gravitational 
collapse and the afterglow. 

If, however, we turn to the detailed temporal structure of the P-GRB and its 
detailed spectral distribution, it is clear that the approximations given in section 
[21 is no longer valid. The detailed description of the formation of the dyadosphere 
as qualitatively expressed in Fig. 1401 is now needed in all mathematical rigour with 
the full development of all its governing equations. Progress in this direction is 
being made at this moment. This situation, however, provides a unique 

opportunity to follow in real time the general relativistic effects of the approach to 
the EMBH horizon as it occurs. In other words all direct general relativistic effects 
of the GRBs are encoded in the fine structure of the P-GRB. For the reasons given 
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Fig. 27. The duration computed with the Tgo criterion is represented as a function of the B 
parameter for three selected EMBH respectively with E^iya = 4.4 X 10^^ erg, E^y^ = 3.1 X 10^'* 
erg, E^ya = 4.1 X 10^* erg going from the lower curve to the upper one. 



in section 1^1 the information on the EMBH mass and charge can only come from 
the short bursts. 

This terminates the injector phase. We now turn to the Beam-Target phase in 
which the ABM pulse collides with the interstellar medium target and the after- 
glow is generated. We shall in the following sections review the basic theoretical 
treatment necessary for the description of these remaining eras and proceed then 
to the confrontation of the EMBH theory with the data. 



12. The era IV: the ultrarelativistic and relativistic regimes in the 
afterglow 

In the introduction we have already expressed the basic assumptions which we have 
adopted for the description of the collision of the ABM pulse with the ISM. In anal- 
ogy and by extension of the results obtained for the PEM and PEMB pulse cases, we 
also assume that the expansion of the ABM pulse through the ISM occurs keeping 
its width constant in the laboratory frame, although the results are quite insensitive 
to this assumption. We assume then that this interaction can be represented by a 
sequence of inelastic collisions of the expanding ABM pulse with a large number of 
thin and cold ISM spherical shells at rest with respect to the central EMBH. Each 
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of these swept up shells of thickness Ar has a mass Ail/jsm and is assumed to be 
located between two radial distances ri and r2 (where r2 — ri = Ar <C ri) in the 
laboratory frame. These collisions create an internal energy AE'int- 

We indicate by Ae the increase in the proper internal energy density due to 
the collision with a single shell and by pB the proper energy density of the swept 
up baryonic matter. This includes the baryonic matter composing the remnant 
around the central EMBH, already swept up in the PEMB pulse formation, and 
the baryonic matter from the ISM swept up by the ABM pulse: 

PB = y • (100) 

Here V is the ABM pulse volume in the comoving frame, Mb is the mass of the 
baryonic remnant and Mism is the ISM mass swept up from the transparency point 
through the r in the laboratory frame: 

Mism = mpnism Y {r^ - rj") , (101) 

where rtip the proton mass and rtism the number density of the ISM in the laboratory 
frame. 

The energy conservation law in the laboratory frame at a generic step of the 
collision process is given by 

PSi7i'Vi + AMismC^ = [pB, + ^^1"°''' + Ae^ yz^Vz. (102) 

where the quantities with the index "1" are calculated before the collision of the 
ABM pulse with an elementary shell of thickness Ar and the quantities with "2" 
after the collision, 7 is the gamma factor and V the volume of the ABM pulse in 
the laboratory frame so that V — 7V. 

The momentum conservation law in the laboratory frame is given by 

Pi3i7it^riVi = ( Pbi — H h Ae I -i^Ur^V^, (103) 

where Ur = — '^ is the radial covariant component of the four- velocity vectonESlSOl 
(see Eq.lEOI)). 



We thus obtain 



. llUryi / AMi^\ 

a - 71 , AMismC^ .-,„^. 



We can use for Ae the following expression 
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because we have assumed a "fully radiative regime" and so -Einti = 0. Substituting 
Ea. (|105|l in Ea. H104|l and applying Ea. H106|) . we obtain: 



Ai^int = Pi3i V'l 1/1 + 271 



PBi Vi \ PBi Vi 



72 = I ■ (108) 



1 + 271^ 



PBl Vl \^ pBj Vl 



These relativistic hydrodynamic (RH) equations have to be numerically integrated. 

These are the actual set of equations we have integrated in the EMBH theory. In 
order to compare and contrast our results with the ones in the current literature, in 
section El ■^^e have introduced the continuous limit of our equations and proceeded 
to have piecewise approximate power law solutions. We examine as well in section 
distill under the above assumptions, the effects of a possible departure from homo- 
geneity in the interstellar medium, still keeping the average density riism — const. 
Although these inhomogeneities are not relevant for the overall behaviour of the 
afterglow which we address here , th ey are indeed important for the actual observed 
flux and its temporal structures^. Also these considerations are affected by the 
angular spreading.EHl 



13. The era V: the approach to the nonrelativistic regimes in the 
afterglow 

The only reason for addressing this last era is that the issue of the approach to 
nonrelativistic behaviour has been extensively discussed in the literature. In our 
treatment these results do not show any particular problems and the relativis- 
tic equations of the previous section continue to hold. In the specific example of 
GRB 991216 we will present in section [T7I some analytic asymptotic expansions of 
these equations. 

This concludes the exposition of the different eras of the EMBH theory. It goes 
without saying that for the description of each era, all the preceding eras must 
necessarily be known in order to determine the space-time grid in the laboratory 
frame and its relation to the arrival times as seen by a distant observer. This is the 
basic message expressed in the RSTT paradigm. 

We can now turn to the comparison of the EMBH theory with the observational 
data. 
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14. The best fit of the EMBH theory to the GRB 991216: the 
global features of the solution 

For reasons already explained in the introduction, we use the GRB 991216 as a 
prototype. We will then later apply the EMBH theory to other GRBs. The relevant 
data of GRB 991216 are reproduced in Fig. |3| the data on the burst as recorded 
by BATSeP and the data on the afterglow from the RXTE satellitJ^ and the 
Chandra satellite j^S see also Halpern et al. (2000).!^ 

The data fitting procedure relies on three basic assumption: 

(1) In the E-APE region, the source luminosity is mainly in the energy band 50- 
300 KeV, so we consider the flux observed by BATSE a good approximation of 
the total flux. 

(2) In the decaying part of the afterglow, we assume that during the R-XTE and 
Chandra observations the source luminosity is mainly in the energy band 2- 
10 KeV, so we can again assume that the flux observed by these satellites is a 
good approximation of the total one. 

(3) We have neglected in this paper the optical and radio emissions, since they are 
always negligible with respect to the X-ray and 7-ray fluxes. In fact, even in 
the latest afterglow phases up to where the X-ray data are available, they are 
one order of magnitude smaller then the X-ray flux. 

These assumptions were initially adopted for the sake of simplicity, but h ave now 
also been justified on the basis of the spectral description of the afterglow.!^ 

As already emphasized in the previous sections, in the EMBH theory there 
are only two free parameters characterising the afterglow: the energy of the dya- 
dosphere, Edya, and the baryonic matter in the remnant of the progenitor star, 
parametrized by the dimensionless parameter B. The location of the remnant has 
been assumed ^ 10^" cm. As discussed in Ruffini et al. (2001a|2 and sectional 
the results are rather insensitive to the actual density a nd lo cation of the baryonic 
component but they are very sensitive to the value of B.EHl 

In Fig. IHl we present the actual first results of fitting our EMBH theory to the 
data from the R-XTE and Chandra satellites, corresponding to selected values of 
Edya and B. There are three distinct features which are clearly evident as a function 
of the arrival time at the detector: an initial rising part in the afterglow luminosity 
which reaches a peak followed by a monotonically decreasing part. 

We have then proceeded to fine tune the two parameters in Fig. |2H1 The main 
conclusions from our model are the following: 

1) The slope of the afterglow in the region where the experimental data are 
present is n = —1.6 and is in perfect agreement with the observational data. The 
index n in this region is rather insensitive to the values of the parameters Edya 
and B. The physical reason for this universality of the slope is rather remarkable 
since it depends on a variety of factors including the ultrarelativistic energy of the 
baryons in the ABM pulse, the assumption of constant average density in the ISM, 
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the "fully radiative" conditions leading predominantly to X-ray emission, as well 
as all the different relativistic effects described in the RSTT paradigm (see also 
section 

2) The afterglow fit does not depend directly on the parameters /Lt, ^ but only 
through their combination Edya- Thus there is a 1-parameter family of values of 
the pair (/i,^) allowed by a given viable value of Edya (see Fig 1161 and sectional). 

3) By fine tuning the parameters of the best fit of the luminosity profile and 
time evolution of the afterglow the following parameters have been found: 

Edya = 4.83 X lO^^erg, B = 3 x 10"^ . (109) 




Fig. 28. Fine tuning of the best fit of the afterglow data of Chandra, RXTE as well as of the 
range of variability of the BATSE data on the major burst by a unique afterglow curve leading to 
the parameter values E^y^ = 4.83 X 10^^erg,B = 3 X 10"'^. 



After fixing in Eq. (|109l) the two free parameters of the EMBH theory, modulo 
the mass-charge relationship which fixes Edya, we can derive all the space-time 
parameters of the GRB 991216 (see Tab. ^ as well as the explicit dependence of 
the gamma factor as a function of the radial coordinate (see Fig. O. 

Of special interest is the fundamental diagram of Fig. ^| Its role is essential 
in interpreting all quantities measured in arrival time (the time of an observer in 
an inertial frame at the detector) and their relations to the ones measured in the 
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laboratory time by an observer in an inertial frame at the GRB source. The two 
times are clearly related by light signals (see Fig. ^ and expressed by the integral 
Eg. 1)37(1 and are also affected by the cosmological expansion (see section 0J). 

15. The explanation of the "long bursts" and the identification of 
the proper gamma ray burst (P-GRB) 

Having determined the two free parameters of the EMBH theory, any other feature 
is a new prediction. An unexpected result soon became apparent, namely that the 
average luminosity of the main burst observed by BATSE can be fit by the afterglow 
curve (see Fig. This led us to the identification of the long bursts observed by 
BATSE with the extended afterglow peak emission (E- APE) . The peak of this E- 
APE occurs at ~ 19.87 s and its intensity and time scale are in excellent agreement 
with the BATSE observations!^ It is clear that this E-APE is not a burst, but is 
seen as such by BATSE due to its high noise threshold.!^ Thus the outstanding 
unsolved problem of explaining the long is radically resolved: the so 

called "long bursts" do not exist, they are just E-APEs (see Fig. I29II . 



80 




1000 



TgQ (seconds) 



Fig. 29. The distribution of the burst durations clearly shows two different classes of events: the 
"short bursts" and the "long bursts" (reproduced from Paciesas et al., 1999*'^). 

We now turn to the most cogent question to be asked: where does one find 
the burst which is emitted when the condition of transparency against Thomson 
scattering is reached? We have referred to this as the proper gam ma ra y burst 
(P-GRB) in order to distinguish it from the global GRB phenomena.^^ We are 
guided in this search by two fundamental diagrams (see Fig. 1121 and Fig. 113(1 : 

(1) In Ruffini, Salmonson, Wilson & Xue (2000 j^it is shown that for a fixed value 
of Edya the value of B uniquely determines the energy Ep-qrb of the P-GRB 
and the kinetic energy Esaryons of the ABM pulse which gives origin to the 
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afterglow (see Fig. I12|l . For the particular values of the parameters given in 
Eq. imi), we find 

Ep-GRB = 7.54 X lif^erg , EBaryons = 9.43 X IQ^'^erg (f f 0) 

and then: 

Ep-GRl^ ^ ^ j-g ^ _ (111) 



EBary 

(2) One important additional piece of information comes from the differences in 
arrival time between the P-GRB and the peak of the E-APE, see Fig.^] Using 
the results of this figure and the numerical values given in Tab.^ we can retrace 
the P-GRB by reading off the time parameters of point 4 in Fig.|21 Transparency 
is reached at 21.57 s in comoving time at a radial coordinate r = 1.94 x 10^'' 
cm in the laboratory frame and at 8.41 x 10~^ s in arrival time at the detector. 

All this, namely the energy predicted in Ea. (|110|l for the intensity of the burst 
and its time of arrival, leads to the unequivocal identification of the P-GRB with the 
apparently inconspicuous initial burst in the BATSE data. We have estimated from 
the BATSE data the ratio of the P-GRB to the E-APE over the noise threshold to 
be ^ 10^^, in excellent agreement with the result in Eq. (|lll|l . see Fig. 1141 

It is important to emphasize that the diagrams in Fig. and Fig. 1121 are not 
universal, but depend on the dyadosphere energy. The corresponding diagrams for 
three selected Edya values {Edya = 5.29 x 10^^ erg, Edya = 4.83 x 10^^ erg and 
Edya — 4.49 X 10^^ erg) are given in Fig. I3l)b where we have plotted the energy of 
the P-GRB and of the E-APE as a function of B. The crossing of the intensity of 
P-GRB and E-APE occurs respectively at Bi = 6.0 x 10-^ B2 = 2.5 x 10"^ and 
B^ = 1.2 X 10~^ where Bi > B2 > B3. In Fig. 130b the same quantities are plotted 
as a function of the baryon mass Mb in units of solar masses and the opposite 
dependence occurs: Mi < M2 < M^. 

The physical reasons beyond these results is the following. We recall that the 
kinetic energy EBaryons and mass Mb of PEMB pulse are 

EBaryor^s = {l-l)MB MB = BEdya (112) 

at the crossing point defined by 

EBaryons ~ Ep-GRB = -^Edya- (H^) 

From these two equations, we obtain 

B^—^-—-:^, (114) 
2(7o - 1) 270 

7o is the Lorentz gamma factor of the PEMB pulse at the transparency point, where 
(see section^} 

/ N ^^B , X 

{n^atr + nB)(JT ^ UbOT = 1, = 1 — — , (115) 
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Fig. 30. a) The same diagram of Fig. |£] is plotted for three different Edya values: Edya = 
5.29 X 10^^ erg (dashed lines), E^y^ = 4.83 X lO'^^ erg (solid Unes) and E^ya = 4.49 X lO'^^ erg 
(dotted lines), b) Same as in a) but plotted as a function of the baryonic mass Mg in units of 
solar masses instead of B. 
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At is the PEMB pulse thickness and Tq the radial position at the transparency 
point. In addition, from the total energy conservation, we have 

(e + ns)7^47rr^A = const., (116) 

where e is the thermal energy of the PEMB pulse. In the regime ^ e, we have 

(117) 

and in the regime ub <C e, we have 

7o~ro. (118) 

Considering the crossing point to occur in the second regime, we obtain at the 
crossing point 

B ^ (EdyayK Mb ~ {Edya)i. (119) 

These results are plotted in Figs.l31b-b. The agreement with the computed results is 
quite satisfactory. The differences can be attributed to the approximation adopted 
in Ea. ljllHI which is modified for high B values. 

The conclusion is that for increasing Edya also the baryonic mass corresponding 
to the cross increases, but in percentage it increases less than Edya- 



16. Considerations on the P-GRB spectrum and the hardness of 
the short bursts 

Regarding the P-GRB spectrum, the initial energy of the electron-positron pairs 
and photons in the dyadosphere for given values of the parameters can be easily 
computed following the work of Preparata, Ruffini & Xue (1998).^^ We obtain 
respectively T — 1.95 MeV and T = 29.4 MeV in the two approximations we have 
used for the average energy density of the dyadosphere (see section IHll. It is then 
possible to follow in the laboratory frame the time evolution of the temperature 
of the electron-positron pairs and photons through the different eras, see Fig. 1321 
The condition of transparency is reached at temperatures in the range of 15 — 55 
KeV at the detector, in agreement with the BATSE results. We emphasize that in 
the limit of B going to 10~* in which the P-GRB coincides with the "short bursts" 
the spectrum of the P-GRB becomes harder in agreement with the observational 
dat#liilH3(seeFig.[7J. 

All the above are average values derived from the two approximations used 
in Fig. El If one wishes to compare the EMBH theoretical results with the fine 
temporal details of the observational data on the P-GRB, a departure from this 
average approach will be needed and the fully time varying relativistic analysis 
outlined in Fig. 021 applies as will be further discussed in section |^ 
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Fig. 31. a) The B values corresponding to the crossings in Fig. lijOb are plotted versus E^y^ (solid 
line). The function B oc E^^^'^ obtained from a qualitative theoretical estimate (see Eg. 11191 ) 
is also plotted (dashed line), b) The Mb values corresponding to the crossings in Fig. l30b are 
plotted versus -Edya (solid line). The function Mg oc i?^^* obtained from a qualitative theoretical 
estimate (see EQ. ill9^ 'l is also plotted (dashed line). 



17. Approximations and power laws in the description of the 
afterglow 

In addition to the BATSE data, there is also clearly perfect agreement with the 
decaying part of the afterglow data from the RXTE and Chandra satellites. 

We can also establish at this point a first set of conclusions on the luminosity 
power law index "n" which is a function depending strongly on the transformation 
t ^ ta (see Fig. I10|l . In the current literature such transformations and the 

corresponding n values are incorrect. Our theoretical value ritheo — —1.6 obtained 
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Fig. 32. The temperature of the pulse in the laboratory frame for the first three eras of Fig. 1 
of Ruffini et al. (2001a)^ is given as a function of the laboratory time. The numbers 1, 2, 3, 4 
represent the beginning and end of each era. The two curves refer to two extreme approximations 
adopted in the description of the dyadosphere. Details are given in Ruffini, Salmonson, Wilson & 
Xue (2000)^0 and in section [HI 



for spherical symmetry for fully radiative conditions and constant density of the 
ISM is in agreement with observed Uobs = —1.616 ± 0.067. No evidence of beaming 
is found in GRB 991216. We shall return to this point in the conclusions. 

An extremely large number of papers in the literature deal with the power law 
index in the afterglow era. This issue has been particularly debated in connec- 
tion with the aim of decreasing the energy requirements of GRBs by the effect of 
beaming. It is currently very popular to infer the existence of bear ning from ths 
direct observations of breakings in the power- law index of the afterglow!^"l^n^^|i^:^l»4|95|06|97|M|08| 
Our aim here is to underline an often neglected point that the power law index of 
the afterglow is the result of a variety of factors including the very different regimes 
in the relation between the laboratory time t and the detector arrival time pre- 
sented in Fig. ^1 No meaningful statements on the values of the power-law index 
of the afterglow can be made having neglected these necessary considerations ex- 
pressed in the RSTT paradigm. This becomes particularly transparent from the 
power law expansion in the semianalytic treatments we present below. It is there- 
fore not so surprising, as we will show in the next session, that the results obtained 
in the EMBH theory differ from the ones in the current literature. 
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17.1. The approximate expression of the hydrodynamic equations 

We proceed to a first approximation and expand Eas. Hl()7lll08|l to second order in 
the quantity 

A A/f. ^2 

<Cl. (120) 



We obtain the following expressions: 

A£;int = (71 - 1) AM«mc' - I nf^\~n\ (^^sm)' , (121) 

Z Mb + jWism 

= - M^l~M + h' nr^l'.f ^2 (AMisn.)' , (122) 

Mb + JWism ^ (A/s + Mi= 

where we set A7 = 72 — 71 and have used the fact that ps^ Vi = {Mb + Mism) (? ■ In 
the limit Ai?int dEint, A7 —> dj, and AMjsm — > c^Msm, neglecting also second 
order terms, where 

dv 

dMism = 'iirr'^mpnisnidr = Anr'^mpnismvdt, v = —, (123) 

and where the ISM number density njsm is assumed for simplicity to be riism = 
Icm^^, we obtain: 

dEi,,t = (7 - 1) dMismc'' , (124) 
dl^- ,/~l dM,sm. (125) 

Mb + JWism 

Eas. (|124lll25|) are limiting cases of Taub's hydrodynamical equations P^EHI^Xhey 
have been at times referred into the GRB literature as the Blandford-McKee equationsP^ 
It is clear that the application of these equations holds if Ea. H12()|l applies. The be- 

P-Bi Vl 



haviour of ^^""i'^ as a function of the radius when Afism ^ Mb is: 



PBi Vl Mb 

The condition Mi^m <C Mb holds for GRB 991216 during the entire evolution of 
the system and so Eq. H120|l is valid (see Fig. I33() . 

Eos. 11241125(1 can be simply solved analytically (see e.g. Blandford & McKee 
1976)Pii«We then have: 

{Mb + M,,^f + C 
{Mb + Mi^nO" - C 

where 



7=777— nr?— ^2 7^' (127) 



C = A/b'^^, (128) 
7o + 1 

where we recall that To and 70 are the radial coordinate and the gamma factor at 
the transparency point and Mb is the initial baryonic mass of the ABM pulse. 
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Fig. 33. The factor is represented as a function of the radial coordinate. It is manifestly 



an increasing function. 



Ea. H127|l is a differential equation for r (t), namely 



cdt 



(Mb + M,,^f + C 



1 -2 



which can be integrated analytically with solutior^^ 



(129) 



Ctq 



2c VC {t-to) = {Mb - m°) (r - ro) 

4 



Cro 



2- 



In 



^1 ° 



B + ^ 



1 



V3 arctan 



B 



(^) 



^ {B + lf 



where m° = -^wmpnismr', 



By/i 

1/3 



arctan 



2- B 

bVs 



(130) 



^ B — *^^o™' ^ and we recall that to is the laboratory 
time at the transparency point. Clearly the fulfilment of Ea. (|120() has to be checked 
to ensure the validity of this solution. 
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17.2. The approximate expression of the emitted flux 



From Eas. H124ll25() . it follows that the emitted flux in the laboratory frame is given 
by (see Fig. EK) 



dE 



inr'^nismmpVj (7 - f ) c^, (131) 
and the corresponding flux in detector arrival time (see Fig. 134b ) by 



47rnismmpC 



dE _\ dt dE 
dti " [dti^ 

7(7-1)^ 



t=t{ti} 



(132) 



For the solution of these equations we distinguish four different phases (A-D). 
The first two correspond to era V. 



Phase A 

Just after the transparency condition is reached, the ISM matter involved is so 
small that we can approximately neglect the Mism term in Ea. H127| and we have: 

7=^7o. (133) 

In the specific case of GRB 991216 we have 70 = 310.1, = 1.94 x lO^'^cm, 
to = 6.48 X 10^ s, ta„ ^ 4.21 X 10"^ s and if^ ~ 8.41 x 10'^ s, where the index 
"o" refers to the quantities at the transparency point. We can then establish the 
following equation describing the ABM pulse motion in this phase: r {t) = vt with 
w ~ c. We can than use the following relation between laboratory time and arrival 
time: 

t = 2jo\. = ^C (134) 

which is in perfect agreement with the full numerical computation (see Fig. I10|l. 
We can substitute these equations into Eas. H131ll33|l . obtaining: 

dE 22 / \ 

— (X7„nismi (135) 
dt 



in laboratory time and 



dE 7o'^ism ud\'^ 



oc 



dti (1 + z 



,3 



tl) (136) 



in arrival time, assuming 7 (7 — 1) ~ 7^. The results of the numerical integration 
of Eas. H104ll05() are in perfect agreement with these approximations (see Fig. I34|l . 
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Fig. 34. a) The GRB flux emitted in laboratory time, b) the flux emitted in the arrival time, 
measured by an observer at rest with respect to the detector (see section ITtI . 



Points P - the two maxima of the energy flux 

Since the contribution of the ISM mass in Eas. (|127l - I128l) can no longer be neglected, 
the value of 7 starts to significantly decrease (see Fig. O and the flux reaches a 
majcimum value. We integrate Ea. (|I3I|) and Ea. (|I33|) using Ea. (|I27l) for 7, assuming 
r {t) = vt with w ~ c and Eq. (|134|l for the relation between the laboratory time and 
the arrival time (see Figs. 1351 llUf) . We can now obtain the point where the emitted 
flux reaches its maximum. In ^neral, the location of the maximum of the flux, 
point P in RufSni et al. (2001a) ,1^ will occur at different events, if considered in the 
arrival time (Pa) or in the laboratory time (Pl)- In this second case, the point Pl 
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is determined by equating to zero the first derivative of Ea. l|131|) . and we have: 



2 Mb 



3 Mi= 



~ 27o, (137) 

Pl 



Me 



620.2. The 

Pl 



which in the case of GRB 991216 gives 7p^ = 206.7 and 

maximum of the observed flux is determined by equating to zero the first derivative 
of Ea. ((nS|) . We obtain: 



5 Me 



IPa -al' 



6 Mi, 



^57o, (138) 

Pa 



Mb 



~ 1550.5. 

Pa 



which in the case of GRB 991216 gives 7p^ ~ 258.4 and 

The results of the numerical integration of Eas. 11041105(1 arc in perfect agree- 
ment with these approximations (see Fig. |^ . 



Phase B - the "golden value" n — —1.6 

In this phase 7 can no longer be considered constant and strongly decreases (see 
Fig. IS)). Mism is increasing, but v is still almost constant, equal to c. As a con- 
sequence, we can still say that r (t) = vt with d = c, but the relation between 
laboratory time and arrival time given in Ea. (|134| ) is no longer valid, and also 
Eq. 1)39(1 is no longer applicable in this phase (see Fig. 110(1 . We can instead write the 
following "effective" relation: 

tocittf', (139) 

which is a result of a best fit of the numerical data in this region. Expanding the 
squares in Eg. 1127(1 . neglecting M?^ with respect to Mg but retaining the terms 
in Mism and assuming 70 S> 1 we obtain: 

where rpj^ and tp^ are the values of r and t at point Pl- Substituting this result 
into Eas. ((131(1 . we obtain the emitted flux in the laboratory frame, given by 

dE 

— « 7p4"ismi^^ , (141) 

and this is in good agreement with the full numerical computation (see Fig. I34(). 

To obtain an analytic formula for the observed flux on the detector, we can still 
try to use the approximate relation between t and given by Ea. l(39() : 

t = 2^it)\ = ^^ti (142) 



where 7 (t) is given by Eg. ((140(1 . We obtain: 
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Using this formula in Ea. H133|) . wc finally obtain: 

''^{<y'^ (144) 



dE -i'pt^ ^ism / ,d\ 
dti °^ (1 + z)- 



(X ni,^{ti) ■ , (145) 



where we again assumed 7 (7 — 1) ~ 7^. This results are not in agreement with the 
observational data, because the power-law index for the observed flux is — 10/7 ~ 
— 1.43, instead of the observed value —1.6. 

This is a confirmation that Ea. H142| l cannot be applied in this phase, as instead 
has been done by many authors in the current literature. We instead have to use 
Ea. (|139|) . In fact, doing so we obtain the correct value: 

dE ^^^-1.6 

The results of the numerical integration of Eas. (|104llU5l) are in perfect agree- 
ment with these approximations (see Fig. 1341) , which implies that the approximate 
Ea. (|124ll25|) can still be used in this regime, but not Ea. (|39|l . which has to be 
replaced by an "effective" local power-law behaviour (see Ea. (|139|) ). 

Phase C 

This new phase begins when 7 has decreased so much that the approximation 
r = ct is no longer valid (see Fig. I35() . In the case of GRB 991216 this happens 
when 7 ~ 3.0, t ~ 1.5 x lO'^ s, ~ 2.9 x 10^ s and r ~ 4.4 x 10^'' cm. In this entire 
phase, r (t) manifests the following behaviour typical of damped motion: 

r{t)^f(l~e-^y (146) 

where f, t* and r are constants that can be determined by the best fit of the 
numerical solution. In the present case of GRB 991216 we obtain: 

f ~ 1.101 X lO^^cm, r ~ 2.072 X lO^s, ~ 4.52 x lOS. (147) 

It is important to note that this interesting behaviour, typical of a damped mo- 
tion, does not lead to any power-law relationship for the emitted flux as a function 
of the laboratory time (see Fig. I34|l . However, if we look at the observed flux as a 
function of the detector arrival time, we see that a power-law relationship still can 
be established, fitting the numerical solution. The result is: 

This quite unexpected result can be explained because the relation between t and 
if depends on r (t) in a nonpower-law behaviour. This fact balances the complex 
behaviour of the emitted flux as a function of the laboratory time, leading finally 
again to a power-law behaviour arrival time. 

In this last phase, however, the fiux decreases markedly, and from the point 
of view of the GRB observations, the most relevant regions are phases A and B 
described above, as well as the peak separating them. 
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Laboratory time (t) (s) 

Fig. 35. The exact numerical solution for r {t) (solid line), together with the line r = ct (dotted 
line) and the fitting function given in EQ. IT461 (dashed hue). 

Phase D 

This last phase starts when the system approaches a Newtonian regime. In the 
case of GRB 991216 this occurs when 7 ~ 1.05, t ~ 5.0 x 10^ s, ~ 3.1 x 10^ s 
and r ~ 1.0 x 10^^ cm. In this phase r it) is again approaching a linear behaviour, 
due to the velocity decreasing less steeply than in Phase C. The emitted flux as a 
function of the laboratory time still does not show a power-law behaviour, while the 
observed flux as a function of detector arrival time does, with an index n = —1.45 
(see Fig. EH). 

18. The power-law index of the afterglow and inferences on 
beaming in GRBs 

The results obtained in the previous sections have emphasized the relevance of the 
proper application of the RSTT paradigm to the determination of the power-law 
index of the afterglow. Particularly interesting is the subtle interplay between the 
different regimes in the relation between the laboratory time and the arrival time 
at the detector clearly expressed by Fig. 1 101 and the corresponding different regimes 
encountered in the first order expansion of the relativistic hydrodynamic equations 
of Taub (1948^^ (see section [TTJ. It is interesting to compare and contrast our 
treatment with selected results of the current literature, in order to illustrate some 
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Table 2. We compare and contrast the results on the power-law index n of the afterglow in the 
EMBH theory with other treatments in the current literature, in the limit of high energy and 
fully radiative conditions. The differences between the values of —10/7 ~ —1.43 (Dermer) and the 
results —1.6 in the EMBH theory can be retraced to the use of the two different approximation 





EMBH theory 


Chiang & Dermer (1999|ifl^ 
Dermer, Chiang & B6ttcher 

B5ttchcr &c Dermer (200oPEl 


Piran /i""nll03l 
Sari & Piran fl999USi] 
Piran (200l|iiJ 


Victri (1997,1"^ 


Halpern et al. (200o(lll 


Ultra-relativistic 


7 = To 
7o = 310.1 

n = 2 
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7 ^ To 

71 ^ 2 






Rclativistic 


^ -3 
3.0 < T < 258.5 
n = -1.6 


„ = = -1.43 


n = = -1.375 




n > -1.47 


Non-rclativistic 


n = -1.36 
1.05 < -y < 3.0 






n = -1.7 




Newtonian 


Ti = -1.45 
1 < -y < 1.05 











relevant points (see Tab.^J. We will consider the results in the literature only with 
reference to the limiting case which we address in our work: the condition of fully 
radiative emission. 

The first line of Tab. [51 describes the ultrarelativistic regime, corresponding to 
an increasing energy flux of the afterglow as a function of the arrival time (phase A 
in previous section). Our treatment and the results in the literature by Dermer et 
al. ^^"^ coincide. They agree as well with the results by Piran et al.^D^^ffiS 

The second line corresponds to the rclativistic regime, in which the energy flux 
of the afterglow, after having reached the maximum (point P in previous section), 
monotonically decreases (phase B in previous section). The dependence we have 
found of the gamma factor on the radial coordinate of the expanding ABM pulse 
does coincide with the one given by Dermer et al. and Piran et al. Our power law 
index n in this regime, which perfectly flts the data, however, is markedly differ- 
ent from the others. Particularly interesting is the difference between our results 
and those of Dermer et al: the two treatments coincide up to the last relation be- 
tween the laboratory time and the arrival time at the detector. As explained in 
Eas. H144ll45|l . the two treatments differ in the approximation adopted in relating 
the laboratory time to the arrival time at the detector, illustrated in Fig. ^| Der- 
mer et al. incorrectly adopted the approximation represented by the lower curve in 
Fig. IIUI and consequently they do not find agreement with the observational data. 
We have not been able to retrace in the treatment by Piran et al. the steps which 
have led to their differen t res ults. Special mention must be made of a result stated 
by Halpern et al. (2000),Elthe last entry in line 2, that an absolute lower limit 
for the power-law index n — 1.47 can be established on theoretical grounds. Such 
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a result, clearly not correct also on the basis of our analysis, has been erroneously 
used ti support the existence of beaming in GRBs, as we will see below. 

The third line in Tab. El is also interesting, treating the nonrelativistic limit 
(Phase C in previous section). This regime has been analysed by Vietri (1997) 
avoiding the exact integration of the equations and relying on simple qualitative 
arguments. These results are not confirmed by the integration of the equations 
we have performed. This is an interesting case to be examined for its pedagogical 
consequences. Having totally neglected the relation between the laboratory time 
and the time of arrival at the detector, which we have illustrated in Fig. EH and 
identifying = t, Vietri reaches a very different power law from our. Moreover, 
his solution brings to an underestimation of the radial coordinate: he estimated a 
radial coordinate of 1.1 x 10^^ cm at = 3.5 x 10"' s, while the exact computation 
shows a result greater than 3.0 x 10^^ cm (see Tab. On the other hand if one 
assumes, from the above mentioned identity = t, t = 3.5 x 10^ s, one obtains a 
gamma factor of ~ 300 (see Tab. ^ in total disagreement with the nonrelativistic 
approximation adopted by Vietri. Quite apart from this pedagogical value, this 
nonrelativistic phase is of little interest from the observational point of view, due 
to the smallness of the flux emitted. 

For completeness, we have also shown our estimates of the index n as the New- 
tonian phase approaches in the last line of Tab. |21 

The perfect agreement between our theoretically predicted value for the power- 
law index, ritheo, and the observed one, Uobs, 

ntheo = -1.6, nobs = -1.616 ± 0.067, (149) 

confirms the validity of our major assumptions: 

(1) The fully radiative regime. 

(2) The constant average density of the ISM {riism = 1 proton/ cm^). 

(3) The spherical symmetry of the emission and the absence of beaming in GRB 991216. 

After the work of Mao & Yi (1994)^^ pointing to the possibility of introducing 
beaming to reduce the energetics of GRBs and after the discovery of the afterglow, 
many articles have appeared trying to obtain theoretical and observation al eviden ce 
for beamed emission in GRBs. The observations have ranged from r adijnnmn to 
opticaf^'^^ o4 110 111 the way to X-rays. Particular attention has been devoted 
to relating the existence of beaming to possible breaks in the light curve slope, 
generally expected at a value of the gamma factor 

where i^o is the beam opening angle. There a re niany articles on this subject; to 
mention only the most popular ones, we recalP^ ^'^ ^ l ^*^ ! ^^ ! ^^ ! Far from having 
reached a standard formulation, these approaches differ from each other in the 
expected time at which the break should take place up to a factor of 20.1^. They 
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differ as well for the opening angle of the beam, up to a factor of Disagreement 
still exists on the number o f breaking points: two in the case ofpa one in the case 
of one again in the case o^^^^^^but differing in position from the one of It 
has also been noticed that other authors have shown thr ough numerical simulations 
that such a transition, if visible at all, is not very sharp!^ 

Ample observational data have been obtained for the GRB 991216, in addition 
to the X-ray band, also in the optical and radio. For the reason mentioned at 
the beginning of section ^1 we only address in this article the problem of the 
7" and the X-ray emission. In that respect, the main arti cle a ddressing the issue 
of beaming in the X-rays for GRB 991216 is the one ofP The key argument 
is based on the theoretical inequality claimed to exist for the power-law index 
n > —1.47 (see above). The fact that the observed X-ray decay rate is found to be 
nobs — 1.616±0.067 is interpreted by the authors as evidence for beaming. Moreover, 
the fact that the decay rate n — —1.6 has been observed before a steepening in the 
optical decay occurred at approximately 1 day of arrival time authorized an even 
more extreme proposal of a narrower beam in the X-rays within the optical beam. 

It is clear from the entire treatment which we have presented and the results 
of the EMBH theory given by ritheo = — 1.6 that there is no evidence for such a 
beaming, as already stated above. The motivation by Halpern et al. (2000) stems 
from the incorrect theoretical assumption of the existence of a lower limit in the 
afterglow power-law index n > —1.47. From our theoretical analysis the existence 
of n = —1.6 is clear proof of isotropic emission in the GRB 991216 and a clear test 
of the complete relativistic treatment of the source. The fact that the break in the 
index should be "achromatic" and the absence of beaming in the X-rays imply an 
absence of beaming also in the optical and radio bands. The observed steepening 
in the optical decay has to find an alternative explanation. Although this is not 
the subject of our present work for the above mentioned reasons, we have found 
interesting the considerations by Panaitescu & Kumar (2001)^^, which find that 
"there are some major difficulties to apply a jet model to GRB 991216". They also 
state, still for GRB 991216, that "the steepening of the optical decay of a few days 
is not due to a jet effect, as suggested by Halpern et al. (2000)!^ but to the passage 
of a spectral break" . 

Concerning our own position on the possibility of beaming in GRBs, we would 
like just to remark that, from a preliminary analysis of beamed emission within the 
EMBH model, we have found some new features which are not encompassed by the 
results in the current literature, and they could become a distinctive signature for 
the discrimination of the existence or nonexistence of beaming. The study of the 
steepening in the optical and radio decay is addressed within the EMBH theory in 
a forthcoming paper.iS^I 
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19. Substructures in the E-APE due to inhomogeneities in the 
Interstellar medium 

The afterglow is emitted as the ABM pulse plows through the interstellar matter 
engulfing new baryonic material. In our previous articles we were interested in 
explaining the overall energetics of the GRB phenomena and in this sense, we have 
adopted the very simplified assumption that the interstellar medium is a constant 
density medium with riism = 1/cm? . Consequently, the afterglow emission obtained 
is very smooth in time. We are now interested in seeing if in this framework we can 
also explain most of the time variability observed by BATSE, all of which except 
for the P-GRB should correspond to the beam-target phase in the IBS paradigm. 

We pursue this treatment still neglecting the angular spreading due to off-axis 
scattering in the radiation of the afterglow. 

Our goal is to focus in this simplified model on the basic energetic parameters 
as well as on the drastic consequences of the space-time variables expressed in the 
RSTT paradigm. 

Having obtained the two results presented in Fig. Inland Fig. 1341 we can proceed 
to attack the specific problem of the time variability observed by BATSE. 

The fundamental point is that in both regimes the flux observed in the arrival 
time is proportional to the interstellar matter density: any inhomogeneity in the in- 
terstellar medium Anism/nism will lead correspondingly to a proportional variation 
in the intensity A/// of the afterglow. This result has been erroneously interpreted 
in the current literature as a burst originating in an unspecified "inner engine" . 

In particular, for the main burst observed by BATSE (see Fig. I36h .^ we have 



There are still a variety of physical circumstances which may lead to such density 
inhomogeneities . 

The additional crucial parameter in understanding the physical nature of such 
inhomogeneities is the time scale of the burst observed by BATSE. Such a burst 
lasts Ata — 20s and shows substructures on a time scale of ~ Is (see Fig. IHfihV In 
order to infer the nature of the structure emitting such a burst we must express 
these times scales in the laboratory time. - Since we are at the peak of the GRB we 
have 7p^ ~ 258.5 (see Ea. (|138|) ') and Ata corresponds in the laboratory time to an 
interval 



which determines the characteristic size of the inhomogeneity creating the burst 
AL - 5.0 X lO^^CTO (see Tab. [2 and Fig.lTUI). 

It is immediately clear from Ea. H151|l and Ea. H152|l that these are the typical 
dimensions and density contrasts corresponding to a small interstellar cloud. As an 
explicit example we have shown in Fig. 1371 the density contrasts and dimensions of 
an interstellar cloud with an average density < n >— 1/cm^. Such a cloud is located 




(151) 



At - 1.0 X lO^s, 



(152) 
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Fig. 36. a) Flux of GRB 991216 observed by BATSE. The enlargement clearly shows the P-GRB.^ 
b) Flux computed in the collision of the ABM pulse with an ISM cloud with the density profile 
given in Fig. 1371 The dashed line indicates the emission from an uniform ISM with n = Icm"^. 
The dotted line indicates the BATSE noise level. 



at a distance of ~ 8.7 x lO^^cm from the EMBH, gives rise to a signal similar to 
the one observed by BATSE (see Fig.Enb)- 

It is now interesting to see the burst that would be emitted, if our present 
approximation would still apply, by the interaction of the ABM pulse with the same 
ISM cloud encountered at later times during the evolution of the afterglow. Fig. 138b 
shows the expected structure of the burst at a distance 4.1 x lO^^cm, corresponding 
to an arrival time delay of ^ 2 days, where the gamma factor is now 7* ~ 3.6. It 
is interesting that the overall intensity would be smaller, the intensity ratio of the 
burst relative to the average emission would remains consistent with Ea. (|151|l . but 
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Fig. 37. The density contrast of the ISM cloud profile introduced in order to fit the observation 
of the burst of GRB991216. The dashed line indicates the average uniform density n = lcm~^. 



the time scales of the burst would be longer by a factor 

shows the corresponding quantities for the same ISM cloud located at a distance 
6.4 X lO^^cm from the EMBH, corresponding to an arrival time delay of ~ 1 month, 
where the gamma factor is ~ 1.5 . 

We return in future worlt^ to examine the angular spreading effects point- 
ing out how they improve the results presented here: the explanation of the time 
variability observed in the so called "long bursts" in the BATSE classification of 
GRBs is confirmed. The smoothness, namely the absence of the above mentioned 
substructures, observed in the latest phases of the afterglow finds as well a most 
natural explanation. 



20. The observation of the iron Unes in GRB 991216: on a possible 
GRB-Supernova time sequence 

We have seen in the previous sections how the time structure of the E-APE gives 
information on the composition of the interstellar matter at distances of the order 
of 5 X 10^^ cm from the source. We would like now to point out that the data on the 
iron lines from the Chandra sa tellite on the GRB 99121^and similar observations 
from other source^^^^^ES make it possible to extend this analysis to a larger 
distance scale, possibly all the way out to a few light years, and consequently probe 
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Fig. 38. (a) Same as Fig. 136b with the ISM cloud located at a distance of 3.17 x lO^^cm from the 
EMBH, the time scale of the burst now extends to ~ 1.58 X IQ^s. (b) Same as (a) with the ISM 
cloud at a distance of 4.71 X Id^"^ cm from the EMBH, the time scale of the burst now extends to 
~ 1.79 X lO^s. 



the distribution of stars in the surroundings of the newly formed EMBH. 

Most importantly, these considerations lead to a new paradigm for the inter- 
pretation of the supernova-GRB correlation.!^ Indeed a correlation between the 

occurrence of GRBs and supernova events exists and has been established by the 
works.116 117 118 119 120 121 103 97 122 

Such an associati on ha s been assumed to indicate that GRBs are generated by 
supernova explosions In turn, such a point of view has implied further conse- 
quences: the optical and radio data of the supernova have been attributed to the 
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GRB afterglow, and many theorists have tried to encompass these data and explain 
them as a genuine component of the GRB scenario. 

We propose instead an alternative point of view implying a very clear distinction 
between the GRB phenomenon and the supernova: if relativistic effects presented in 
the RSTT paradigm are properly taken into account, then a kinematically viable ex- 
planation can be given of the supernova-GRB association. We still use GRB 991216 
as a prototypical case. 

The GRB-Supernova Time Sequence paradigm, which we have indicated for 
short as GSTS paradigm,!^ states that: A massive GRB-progenitor star Pi of mass 
Ml undergoes gravitational collapse to an EMBH. During this process a dyadosphere 
is formed and subsequently the P-GRB and the E-APE are generated in sequence. 
They propagate and impact, with their photon and neutrino components, on a second 
supernova-progenitor star P2 of mass M2. Assuming that both stars were generated 
approximately at the same time, we expect to have M2 < Mi. Under some special 
conditions of the thermonuclear evolution of the supernova-progenitor star P2, the 
collision of the P-GRB and the E-APE with the star P2 can induce its supernova 
explosion. 

Especially relevant to our paradigm are the following data from the Chandra 
satelliter^S 

(1) At the arrival time of 37 hr after the initial burst there is evidence of iron 
emission lines for GRB 991216. 

(2) The emission lines are present during the entire observation period of 10'* s. The 
iron lines could also have been produced earlier, before Chandra was observing. 
Thus the times used in these calculations are not unique: they do serve to 
provide an example of the scenario. 

(3) The emission lines appear to have a peak at an energy of 3.49 ± 0.06 keV 
which, at a redshift z = 1.00 ± 0.02 corresponds to an hydrogen-like iron line 
at 6.97 keV at rest. This source does not appear to have any significant motion 
departing from the cosmological flow. The iron lines have a width of 0.23 keV 
consistent with a radial velocity field of 0.1c. The iron lines are only a small 
fraction of the observed flux. 

On the basis of the explicit computations of the different eras presented in the 
above sections, we make three key points: 

(1) An arrival time of 37 hr in the detector frame corresponds to a radial distance 
from the EMBH travelled by the ABM pulse of 3.94 x 10^''' cm in the laboratory 
frame (see Tab.^. 

(2) It is likely that a few stars are present within that radius as members of a cluster. 
It has become evident from observations of dense clust ers o f star- forming regions 
that a stellar average density of typically lO^pc^'^ [i2M should be expected. 
There is also the distinct possibility for this case and other systems that the 
stars Pi and P2 are members of a binary system. 



On the structure of the burst and afterglow of Gamma-Ray Bursts I: the radial approximation 81 



(3) The possible observations at different wavelengtlis of tlie supernova crucially 
depend on the relative intensities between the GRB and the supernova as well 
as on the value of the distance and the redshift of the source. In the present case 
of GRB 991216, the expected optical and radio emission from the supernova 
are many orders of magnitude smaller than the GRB intensity. The opposite 
situation will be encountered in GRB 980425.^3 

In order to reach an intuitive understanding of these complex computations we 
present a schematic very simplified diagram (not to scale) in Fig. 1391 

We now describe the sequence of events and the specific data corresponding to 
the GSTS paradigm: 

(1) The two stars Pi and P2 are separated by a distance Dp.^ = 3.94x 10^^ cm in the 
laboratory frame, see Fig. Both stars are at rest in the inertial laboratory 
frame. At laboratory time t — Q and at comoving time t = 0, the gravitational 
collapse of the GRB-progenitor star Pi occurs, and the initial emission of grav- 
itational radiation or a neutrino burst from the event then synchronizes this 
event with the arrival times = at the supernova-progenitor star P2 and 

= for the distant observer at rest with the detector. The electromagnetic 
radiation emitted by the gravitational collapse process is in stea d practically 
zero, due to the optical thickness of the material at this stage,!^ see Tab.Q] 

(2) From Tab. ^ at laboratory time ti = 6.48 x 10^ s and at a distance from 
the EMBH of Di ~ 1.94 x 10^'* cm, the condition of transparency for the 
PEMB pulse is reached and the P-GRB is emitted (see section [HI. This time is 
recorded in arrival time at the detector t'^_^ = 8.41 x 10~^ s, and, at P2, at ta^ = 
4.20 X 10"^ s. The fact that the PEMB pulse in an arrival time of 8.41 x 10"^ s 
covers a distance of 1.94 x 10^^ cm gives rise to an apparent "superluminal" 
effect. This apparent paradox can be straightforwardly explained by introducing 
an "effective" gamma factor.^ 

(3) At laboratory time t — 1.73 x 10^ s and at a distance from the EMBH of 
5.18 X 10^^ cm in the laboratory frame, the peak of the E-APE is reached 
which is recorded at the arrival time ta = 9.93 s at P2 and = 19.87 s at the 
detector. This also gives rise to an apparent "superluminal" effect. 

(4) At a distance Dp.^ = 3.94 x 10^^ cm, the two bursts described in the above 
points 2) and 3) collide with the supernova-progenitor star P2 at arrival times 
ta^ — 4.20 X 10~^s and ta = 9.93 s respectively. They can then induce the 
supernova explosion of the massive star P2 ■ 

(5) The associated supernova shell expands with velocity 0.1c. 

(6) The expanding supernova shell is reached by the ABM pulse generating the 
afterglow with a delay of = 18.5 hr in arrival time following the arrival 
of the P-GRB and the E-APE. This time delay coincides with the interval 
of laboratory time separating the two events, since the P2 is at rest in the 
inertial laboratory frame.^ The ABM pulse has travelled in the laboratory 
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Fig. 39. A qualitative simplified space-time diagram (in arbitrary units) illustrating the GSTS 
paradigm. The EMBH, originating from the gravitational collapse of a massive GRB-progenitor 
star Pi , and the massive supernova-progenitor star P2-neutron star (P2-NS) system, separated by 
a radial distance Dp^ , are assumed to be at rest in in the laboratory frame. Their worldlines are 
represented by two parallel vertical lines. The supernova shell moving at 0.1c generated by the 
P2-NS transition is represented by the dotted line cone. The solid line represents the motion of the 
pulse, as if it would move with an "effective" constant gamma factor 71 during the eras reaching the 
condition of transparency. Similarly, another "effective" constant gamma factor 72 < 71 applies 
during era IV up to the collision with the P2-NS system. A third "effective" constant gamma 
factor 73 < 72 occurs during era V after the collision as the nonrelativistic regime of expansion is 
reached. The dashed lines at 45 degrees represent signals propagating at speed of light. 

frame a distance Dp^ — _Di ~ Dp^ = 3.94 x 10^^ cm in a laboratory time 
f2 — — f2 = 1-32 X 10'' s (neglecting the supernova expansion). 

The collision of the pulse with the supernova shell occurs at 7 ~ 4.0. By this time 
the supernova shell has reached a dimension of 1.997 x lO""^^ cm, which is consistent 
with the observations from the Chandra satellite. 
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In these considerations on GRB 991216 the supernova remnant has been as- 
sumed to be close to but not exactly along the line of sight extending from the 
EMBH to the distant observer. If such an alignment should exist for other GRBs, 
it would lead to an observation of iron absorption lines as well as to an increase 
in the radiation observed in the afterglow corresponding to the crossing of the su- 
pernova shell by the ABM pulse. In fact, as the ABM pulse engulfs the baryonic 
matter of the remnant, above and beyond the normal interstellar medium baryonic 
matter, the conservation of energy and momentum implies that a larger amount 
of internal energy is available and radiated in the process (see section IT^ . This 
increased energy-momentum loss will generally affect the slope of the afterglow de- 
cay, approaching more rapidly a nonrelativistic expansion phase (details are given 
in section [T7|l . 

It is quite clear that as soon as the relativistic transformations of the RSTT 
paradigm are duly taken into account, the sequence of events between the super- 
nova and the GRB occurrences are exactly the opposite of the one postulated in the 
so-called "supranova" scenario. I^^^^^l This can be considered a very appropri- 
ate pedagogical example of how classical nonrelativistic applied to ultrarelativistic 
regimes can indeed subvert the very causal relation between events. 

If we now turn to the possibility of dynamically implementing the scenario, there 
are at least three different possibilities: 

(1) Particularly attractive is the possibility that a massive st ar has rapidly 
evolved during its thermonuclear evolution to a white dwarf.^^ It it then suf- 
ficient that the P-GRB and the E-APE implode the star sufficiently as to reach 
a central density above the critical density for the ignition of thermonuclear 
burning. Consequently, the explosion of the star P2 occurs, and a significant 
fraction of a solar mass of iron is generated. These configurations are currently 
gene rally considered precursors of some type I supernovae (see e.g. Filippenko, 
1997j I and references therein). 

(2) Alternatively, the massive star P2 can have evolved to the condition of being 
close to the point of gravitational collapse, having developed the formation of an 
iron-silicon core, type II supernovae. The above transfer of energy momentum 
from the P-GRB and the E-APE may enhance the capture of the electrons on 
the iron nuclei and consequently decrease the Fermi energy of the core, leading 
to the onset of gravitational instability (see e.g. Bethe, 1991,"'^^''' p. 270 and 
followings). Since the time for the final evolution of a massive star with an 
iron-silicon core is short, this event requires a well tuned coincidence. 

(3) The pressure wave may trigger massive and instantaneous nuclear burning pro- 
cess, with corresponding changes in the chemical composition of the star, lead- 
ing to the collapse. 

The GSTS paradigm has been apphed to the case of the GRB 980425 - SN1998bw 
which, with a red shift of 0.0083, is one of the closest and weaker GRBs observed. 
In this case, the radio and the optical emission of the supernova is distinctively 
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observed. For this particular case, the EMBH appears to have a significantly lower 
value of the parameter ^ and the validity of the GSTS paradigm presented here is 
confirmed!"^ 



21. General considerations on the EMBH formation 

Before concluding let us consider the problem of the EMBH formation. Such a 
problem has been debated for many years since the earliest discussions in 1970 in 
Princeton and has been finally cla rifie d and addressed in general terms to justify 
the plausibility of the hypothesis in.^^ There has been a basic change of paradigm. 
All the considerations on the electric charge of st ars w ere traditionally directed, 
following the classical work by Shvartsman (197oj^Sl all the way to the funda- 
mental book by Punsly,^^ to the presence of a net charge on the star surface in 
a steady state condition. The star can be endowed of rotation and magne tic fi eld 
and surrounded by plasma, hke in the case of Goldreich & JuHan (1969) or, 
in the case of absence of both magnetic field and rotation, the electrostatic pro- 
cesses can be relat ed to the depth of the gravitational well, like in the treatment of 
Shvartsman (1970).II^ However, in neither cases it is possible to reach the condi- 
tion of the overcritical field needed for pair creation nor have the condition of no 
baryonic contamination discussed in sections [3 El and essential for the dyadosphere 
formation. The basic conceptual point is that GRBs are maybe the most violent 
transient phenomenon occurring in the universe and so the condition for the dya- 
dosphere creation have to be searched in a transient phenomenon. The solution is 
related to the most transient phenomenon occurring in the life of a star: the process 
of gravitational collapse. 

Having acquired such a fundamental understanding, the next step is to estimate 
the amount of polarization needed in order to reach the fully relativistic condition 

= 1 . (153) 

Recalling that the charge to mass ratio of a proton is qp/ {rnp\/(Pj = 1.1 x 10"'^*, 

it is enough to have an excess of one quantum of charge every 10^* nucleons in 
the core of the collapsing star to obtain an extreme EMBH after the occurrence of 
the gravitational collapse. Physically this means that we are dealing with a process 
of charge segregation between the core and the outer part of the star which has 
the opposite sign of net charge in order to enforce the overall charge neutrality 
condition. We here emphasize the name "charge segregation" instead of the name 
"charge separation" in order to contrast a very mild charge surplus created in dif- 
ferent part of the star, keeping the overall charge neutrality, from the much more 
extreme condition of charge separation in which all the charges of the atomic com- 
ponent of the star are separated. It is indeed reassuring that such a core, endowed 
with charge segregation, is indeed stable with respect to the Fermi-Chandrasekhar 
criteria for the stability of self-gravitating stars duly extended from the magnetic 
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to the electric case: the electric energy of such a core is consistently smaller than 
its gravitational energy!^ 

Such a condition of charge segregation between the core and the oppositely 
charged star surface layer can be reached under a very large number of physical 
conditions. We consider, for simplicity, one of the oldest example: the one of a star 
endowed with both a magnetic field and rotation. It is proved that a typical mag- 
netic field expected for the ISM is Bo ~ 10" 

-5 cjlMiWe further assume, consistently 
with the data which we have acquired and verified in the present article (see sec- 
tions El EIj that also in the galaxy where GRB 991216 occurred the ISM has an 
average density oiuism = \ proton / cm^ . From this value of density we have that an 
ISM cloud with mass M ^ IQMq occupies a sphere of radius Ro ~ 1.4 x 10^^ cm. 
If this sphere collapse to a star with radius R = Rq , from the flux conservation we 
obtain that it is enough for this star to rotate with the most reasonable angular 
speed 



" ^ RqR^oBo ^^^^^ 

to conclude that the progenitor star core is endowed of a charge to mass ratio equal 
to In the extreme case of Ea. (|153ll we have ^ = 1 and so the angular speed is 
fl ~ 1.1 X 10^^ rad/s — i. e. o ne round in 1.5 /ir — and correspondingly we have 
smaller il values for ^ < iP^ Clearly the overall neutrality is guaranteed by the 
oppositely charged baryonic matter which is the one measured by the B parameter 
in the EMBH model (see sections CHHJ- The smallness of the B value clearly points 
to the absence of an extended envelope of the progenitor star. 

The formation process of such an electromagnetised progenitor star will be 
clearly affected by the presence of differential rotation, the consequent amplifica- 
tion of the magnetic field and a variety of magnetohydrodynamical problems which 
will affect somewhat the simplicity of the heuristic Ea. H154|l . Similarly the process 
of gravitational collapse of such a progenitor star endowed with rotation will lead 
to complex phenomena of "gravitationally induced electrornagnetic radiation'^^21 
and of "electromagnetically induced gravitational radiation'^-"^ which will tend to 
reduce both the eccentricity and the angular velocity of the collapsing core. The 
general outcome of gravitational collapse will be a Kerr-Newmann spacetime. It is 
interesting that such a general case will break the degeneracy in (/i, ^) described in 
section 1^1 -'In this article we have addressed the much simpler case of a solution in 
which (cL) I {GM'^^ <^ 1 and the treatment can be well approximated by a collapse 
described by a Reissner-Nordstrom geometry. 

In addition to this scenario, based on the role of magnetic field and rotation, we 
are as well pursuing the possible generation of the charge segregation by quantum 
effects at the surface of the almost Fermi degenerate core. This most straightforward 
analysis also leads to a Reissner-Nordstrom geometry. 

In both these cases the Reissner-Nordstrom geometry appears indeed to be the 
relevant model for GRB 991216 as discussed in the previous sections. We shall return 
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to non spherical configuration in forthcoming publications and/or when requested 
by observational evidence @1 

Turning now from this general scenario to a more detailed analysis of a Reissner- 
Nordstrom geometry, some preliminar y ne cessary steps have to be accomplished. In 
Cherubini, Ruffini & Vitagliano 

(2002 jZ2 

we have considered the gravitational col- 
lapse of a charged spherical shell with selected boundary conditions: either starting 
from infinite distance with a null or non null kinetic energy, or imploding from a 
finite distance initially at rest. A new analytic solution has been obtained for such 
a boundary condition, corresponding both to a collapse into an already formed 
EMBH or to a collapse in Minkowsky space. In both cases we have followed the 
process of gravitational collapse all the way to the self closure of the shell by the 
formation of an horizon. 

Using this analytic solution it has been possible to clarify the independent phys- 
ical components, contributing to the formation of the EMBH irreducible mass.L^ 
Surprisingly, the irreducible mass does not directly depend on the electromagnetic 
energy of the imploding shell: it is uniquely a function of the initial baryonic mass, 
of its gravitational energy and of the kinetic energy of the implosion. The electro- 
magnetic energy is stored around the EMBH and can be extracted by two very 
different process as a function of the electromagnetic field strength, a) When the 
electric field on the collapsing shell is smalle r than £r , the process of energy extrac- 
tion occurs in the effective EMBH ergospher a sequence of discrete high 
energy events, with energy up to 10^^-10^^ eV. Such s ource s can be of relevance 
for the explanation of the ultra high energy cosmic rays.EEl b) When the electric 
field on the collapsing shell is larger than £c, the conditions relevant to the present 
article are fulfilled. The energy extraction process occurs in the dyadosphere and a 
much larger number of electron and positron pairs are created with typical energies 
of the order of 10 MeV which are relevant for the process considered in the present 
paper. 

It is interesting that the clarification obtained in Ruffini & Vitaghano (2002aj^ 
has allowed a deeper understanding of the essential role of the gravitational and 
kinetic implosion energies and the storage of the electromagnetic energy in the entire 
region surrounding the EMBH horizon. It has been shown in Ruffini & Vitagliano 
(2GG2b)'^4 ti^at the central point can simply be summarized: the Coulomb repulsion 
of the collapsing matter reduces the kinetic energy of implosion leading to a smaller 
value of the irreducible mass and consequently to a larger value of the extractable 
energy. 

Having so established and clarified the basic conceptual processes of the ener- 
getic of the EMBH, we are now ready to approach, using the new analytic solu- 
tion obtained, the dynamical process of vacuum polarization occurring during the 
formation of an EMBH as qualitatively represented in Fig. QHI The study of the 
dyadosphere dynamical formation as well as of the electron-positron plasma dy- 
namical evolution will lead to the first possibility of directly observing the general 
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relativistic effects approaching the EMBH horizon. 

Before closing we would like to emphasize once more a basic point: all the 
considerations presented in the description of the preceding eras are based on the 
approximations in the description of the dyadosphere presented in section El This 
treatment is very appropriate in estimating the general dependence of the energy 
of the P-GRB, the kinetic energy of the ABM pulse and consequently the intensity 
of the afterglow, as well as the overall time structure of the GRB and especially the 
time of the release of the P-GRB in respect to the moment of gravitational collapse 
and its relative intensity with respect to the afterglow. If, however, is addressed 
the issue of the detailed temporal structure of the P-GRB and its detailed spectral 
distributio n, t he above dynamical considerations on the dyadosphere formation 
are needed.L^ In turn, this detailed analysis is needed if the general relativistic 
effects close to the horizon formation have to be followed. As expressed already 
in section. 1111 all general relativistic quantum field theory effects are encoded in 
the fine structure of the P-GRB. As emphasized in section O the only way to 
differentiate between solutions with same Edya but different EMBH mass and charge 
is to observe the P-GRBs in the limit S ^ 0, namely, to observe the short GRBs. 



Dyadosphere 




space 



Fig. 40. Space-time diagram of the collapse process leading to the formation of the dyadosphere. 
As the collapsing core crosses the dyadosphere radius the pair creation process starts, and the pairs 
thermalize in a neutral plasma configuration. Then also the horizon is crossed and the singularity 
is formed. 
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22. Conclusions 

This paper is a consequence of a process of revision, r eanalYsis arid fu rther expansion 
of all the results presented in the previous article d I I * * I I In this process 
all the numerical estimates have been reanalysed for consistency. Most noticeable 
we have found a missing factor of 27r in the values of the experimental data on 
GRB 991216 we had used in our previous works which led to an underestimate of 
the total energy of the source. 

This process of revision, far from being just a detailed computational verifica- 
tion, has given us the opportunity to rethink the entire GRB process in an unitary 
description starting from the moment of gravitational collapse all the way up to 
the latest phases of the afterglow and in identifying the three fundamental actors 
of the GRB phenomenon: 

(1) Edya- Having reanalysed in section |2 the physics of the dyadosphere we have 
pointed out in Fig. ^] that the same value of Edya can be obtained from an 
entire family of {^J.^C) parameters (i.e. Edya is degenerate in (/x,^)). We have 
then shown in the reexamination of all the GRB eras that all the results de- 
pend only on the value of Edya and not on the particular value of (/i, ^) (see 
sections I7I81 121131) . The only exception to this occurs in the era I (see section|Bl) 
which is the only one relevant for short GRBs. 

(2) B. The crucial role played by the baryonic remnant of the progenitor star in 
determining the relative intensity ratio and the time delay between the P-GRB 
and the E-APE has been summarized already in the two Figs. I12MT^ in the 
introduction. 

(3) ISM. The density n^sm of the interstellar medium and its inhomogeneities ap- 
pears to have a fundamental role in the intensity and the temporal substruc- 
tures of the E-APE and the afterglow. In order to identify such a crucial role, 
however, the correct relativistic space-time relations expressed by the RSTT 
paradigm are needed as amply exemplified in sections 1 17H 191 

The observational data agree with the predictions of the model on: 

1) the intensity ratio, 1.58 x 10~^, between the P-GRB and the E-APE, which 
strongly depends on the parameter B, 

2) the absolute intensities for both the P-GRB and the E-APE, respectively 7.54 x 
10^1 erg and 4.75 x lO^^^ 

3) the arrival time of the P-GRB and the peak of the E-APE, respectively 8.41 x 10~^ 
s and 19.87 s. 

These results can certainly be considered the greatest success of the EMBH 
theory. 

Before closing, we like to draw some specific conclusions based on the first 
fundamental parameter of the EMBH theory: Edya- It is clear that Edya is the 
fundamental parameter which determine the general energetic requirements of the 
GRB 991216. This energetics strongly depends on the possible existence or absence 
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of beaming in the radiation process. In turn, as expressed in sections 117m 81 the 
presence of beaming is led back to the power-law index n of the afterglow. The 
general conclusions reached on Edya can be summarized as follow: 

1.1 The value of n is a function of the transformation between t ^ ta ^ (see 
section I17II . These transformations are a function of the entire relativistic regime 
of the world line of the source (see section Q . By systematically neglecting this 
information the current works in the GRBs literature have obtained incorrect n 
values (see Tab.|21and section IT^ . 

1.2 The value ntheo = —1.6, which we have obtained within the EMBH theory, in 
the region of interest for the observations, based on the assumptions of spherical 
symmetry, fully radiative condition in the emission process and constant density of 
ISM, is in agreement with observed Uobs = —1.616 ±0.067 (see sections [T7lll8|l . No 
evidence of beaming is therefore found is GRB 991216 (see section IT^ . 

1.3 For GRB 991216 E^ya = 4.83 x 10^^ gj.g jg fo^j^j the EMBH theory. This 
value is systematically lar ger t han the ones quoted in th e cu rrent literature by 
Panaitescu & Kumar (2001^^ and by Halpern et al. (2000 due to the fact that 
they respectively consider beaming angles of 3° — 4° and 6°. These considerations 
have been shown to be untenable in sectional There is still a difference of ~ 28% 
between the total energy implied by the EMBH theory (4.83 x 10^^ erg) and the 
value quoted by Halpern {E^ya = 6.7 x 10^'^ erg) in the case of spherical emission. 
We trust that this is a consequence of the underlying assumption of the spectral 
distribu tion of the radiation assumed by Halpern et al. (2000)^ (see e.g. Frail et 
al., 2001^23), which should be reassessed on the ground of our theoretical results. 
See also paper lI.El 

We turn now to the second fundamental parameter of the EMBH theory: B. 
It is essential in explaining the difference between the so called "long bursts" and 
"short bursts" (see sections ^2 115|l . The general conclusions reached on B can be 
summarized as follows: 

2.1 The most general GRB contains three different components: the P-GRB, the 
E-APE and the rest of the afterglow. The ratio between the P-GRB and the E- 
APE intensity and their temporal separation is a function of the B parameter (see 
Figs. nHIT^ . The best fit is obtained for B = 3.0 x lO^^ (ggg sectionlHI). We recall 
that in the present case for B < 2.5 x 10^^ the energy of the P-GRB would be larger 
than the one of the E-APE and the energy of the dyadosphere would be mainly 
emitted in what have been called the "short bursts", while for B > 2.5 x 10^^ the 
energy of the E-APE would predominate and the energy of the dyadosphere would 
be mainly carried by the ABM pulse and emitted in the afterglow. 
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2.2 The difficulties encountered by all theoretical models, through the years, in 
order to explain the so called "long bursts" are resolved in a drastic way (see 
section EJ- The so called "long bursts" are not bursts at all. They represent just 
the E-APE which was interpreted as a burst only due to the noise threshold in the 
BATSE observations (see Fig. EJ- The E-APE is emitted at distances from the 
EMBH in the range 1.0 x 10^^ ^ 1.0 x 10^^ cm, see Tab. Q namely well outside 
the size of the progenitor star and already deep in interstellar space. The fact 
that the crossing of such distance, which is a typical dimension of an interstellar 
cloud, appears to occur in arrival time in only ~ 100 seconds is perfectly explained 
by the relativistic transformations encoded in the RSTT paradigm corresponding 
to a gamma factor between 100 and 300 (see section 0] and Tab. QJ. This effect 
would be interpreted within a classical and incorrect astronomical picture by a 
"supcrluminal" behaviour propagating at ~ 3.6 x lO^'c (see Tab.^. 



2.3 In the limit B ^ the entire dyadosphere energy is emitted in the P-GRB. 
These events represents the "short bursts" class, for which the afterglow intensity 
is smaller than the P-GRB emission and below the actual observational limits (see 
section [TT|l . It is interesting that the proposed differentiation between the "short 
bursts" and "long bursts" within the EMBH theory is merely due to the amount 
of baryonic matter in the remnant, described by the B parameter, and totally 
independent from the process of gravitational collapse which is clearly identical in 
both cases. This explains at once the recently found conclusion that the distribution 
of short and long GRBs have essentially the same characteristic peak luminosity.l^ 
Also the result expressed in Fig. |3 that the average temperature corresponding to 
the P-GRB emission does increase for decreasing values of the B parameter can 
explain the observed fact that the "short bursts" , which are obtained in the limit 
B ^ 0, are systematically harder than "long bursts" I^Sl 

Finally, the EMBH theory offers an unprecedented tool in order to map with 
great accuracy all the matter distribution around the newly formed EMBH from 
the horizon all the way to the ISM. This concept was pioneered by Dermer & 
Mitman (1999^P^ who proposed to use GRB sources as "tomographic images of 
the density distributions of the medium surrounding the sources of GRBs" . It is 
important to emphasize that the very precise reading of the matter distribution 
encoded in the data of the P-GRB, the E-APE and the afterglow in GRB 991216 is 
in mark ed disagre ement with the matter distribution postulated by the "coUapsar" 
scenario .1^3^^1221 This conclusion is evidenced not only by the absence of beaming 
already mentioned above, but also for the paucity of the baryonic matter encoun- 
tered by the PEM pulse in its way out from the EMBH. There is no evidence for 
the presence either of a baryonic disk component nor of a conspicuous baryonic 
remnant. We actually have B — 3.0 x 10^^. The general conclusions reached on 
this topic can then be summarized as follow: 
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3.1 Starting from the inside out we have, from the electrodynamics of the dyado- 
sphere, that the baryonic contamination in that region has to be much smaller than 
10®g/cm'^ (see section EJ. This condition can be achieved in the formation of an 
EMBH. The same electrodynamical process would be hampered in the formation 
of a neutron star due to the high density and baryonic contamination. Among all 
process of discharge of the overcritical electromagnetic field the pair creation is the 
most effective one due to the very short time scale of the order of — ~ 10^^^ 
seconds (see section EJ. 

3.2 Unlike the case of formation of a neutron star, the mass of the remnant of the 
progenitor star is very small indeed. This mass, determined by -B, is very accurately 
inferable from the relative intensity and temporal distance between the P-GRB and 
the E-APE (see above). In the present case we have Mb ~ 8.1 x 10"'' Mq. The pres- 
ence of the remnant is also important for guaranteeing the overall charge neutrality 
of the system formed by the oppositely charged collapsing core and the remnant. It 
has been pointed out in section [?Tl that this condition of charge separation between 
the collapsing core and the remnant occurs only during the relevant part of the 
gravitational collapse process which, we recall, for a IOA/0 is of the order of 30 
seconds. 

3.3 The structure of the E-APE and the afterglow gives as well an unprecedented 
tool in order to estimate the average density and filamentary structure in the ISM: 
the structures down to a fraction of seconds observed in the E-APEs, the so called 
"long bursts" of the current literature (see Fig. OJ, can be used in order to map 
the filamentary structure as well as the size of interstellar clouds surrounding the 
EMBH (see section ITi^ . When all the geometrical and relativistic effects are duly 
taken into account the intensity and the average profile of the E-APE and of the 
afterglow point to an average value of the ISM density riism ~ Iproton / cm^ in very 
good agreement with a large variety of independent estimates. The very late phases 
of the afterglow gives information of the induced supernova collapse (see section EHl 
and Ruffini et al., 2001c^) which will be addressed in a forthcoming publication. 
Since now we can assert that the correct space-time se quence ba sed on the RSTT 
paradigm is in contrast with the "supranova" scenariJ^^^^^^ which was based 
on a nonrelativistic consideration in ultrarelativistic regimes (see section 1201) • 

This concludes the presentation of the basic model which is now ready to be 
applied to additional sources. If we look to the future we can see three main topics 
to be addressed with special attention: 

(1) We have performed a more detailed description of beaming, of the angular 
spreading and of the spectral properties which is going to be the subject of 
paper II.I^ Since now, we can assert that this more detailed treatment supports 
all general conclusions obtained in the present paper. 

(2) If one is interested in the detailed effects of general relativity and relativistic 
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field theory, all the attention should be directed to the structure of the short 
bursts. This needs the development of detailed theoretical works on the ap- 
proach of the horizon of the black hole and the associated electrodynamical 
process. From the description presented in this paper of an already formed and 
averaged dyadosphere (see Fig. 117(1 we have to move to the treatment of its 
dynamical formation (see Fig. 140(1 . Such an analysis describing the approach to 
the formation of the horizon of t he EMB H. within the EMBH theory, is in ad- 
vanced phase of development*^^ 7 3174 1 7^?! gee Fig. jlHl Some preliminary results 
have appeared in Bianco, Rufiini & Xue (2001)!^ 
(3) From the observational point of view, the detailed observations of the yet un- 
explored region in the range up to 10^ seconds in Fig. 1141 and the corresponding 
observations of the "short bursts" by a new class of space missions with higher 
sensitivity than the BATSE instrument appear to be of great importance. Such 
observations should allow to directly observe for the first time the general rel- 
ativistic and extreme quantum field theory effects connected to the process of 
formation of the EMBH. It can be of some interest to explore the possibil- 
ity of obs ervin g in these regimes the "gravitationally induced electromagnetic 
radiation'^E^and the "electromagnetically induced gravitational radiation''»i23l 
phenomena. 
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